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ABSTRACT 

Our  extensions  to  conical  orlficta  of  Clausing  'a  analysis  of  angular 
distribution  of  molecules  effusing  from  cylindrical  orifices  has  resulted 
in  numerical  values  for  transmission  probabilities  and  recoil-force  correction 
factors  which  are  tabulated.  With  these  results,  it  is  demonstrated  that 
the  opLiiuuiu  orifice  geometry  for  ( l)  recoil  force  measurements  is  a  diverging 
conical  orifice  with  semi-apex  angle  of  30°,  (.?)  delivery  of  maximum  fraction 
of  effusing  molecules  onto  (or  into)  a  target  (aperture)  is  a  long  cylindrical 
orifice. 

Modifications  to  the  angular  distribution  apparatus  are  described  and 
experimental  results  given  for  four  orifices  and  two  gases  over  the  pressure 
range  5  to  900  dyn/cm3.  The  most  interesting  aspect  ol  the  results,  one 
apparently  not  previously  noticed  in  angular  distribution  results,  is  the 
presence  of  maxima  and  minima  in  plots  of  A  vs.  9,  where  A  is  the  (experimental 
value  -  theoretical  value)  of  Pq,  the  fraction  of  effusing  molecules  which 
flow  per  steradian  at  angle  G  from  the  orifice  axis.  These  maxima  and  minima 
have  been  correlated  with  the  relauive  contribution  from  the  orifice  wall 
to  the  total  flux  at  angle  G. 

Additional  results  for  experimental  transmission  probabilities  of  orifices 
determined  by  the  Multicell  technique  are  generally  in  agreement  with  theo¬ 
retical  values  within  2  to  3$. 

The  Miker  technique  for  simultaneous  determination  oi  vapor  pressure  by 
rate  of  effusion  and  by  recoil  force  measurements  has  been  refined  to  the  point 
that  recoil  force  data  are  as  reproducible  as  rate  of  effusion  measurements. 
Several  sources  of  spurious  recoil  force  have  been  identified  and  eliminated. 

A  new  furnace  and  a  modified  automatic  control  system  for  the  microbalance 
are  described. 
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INTRODUCTION 

The  Knudsen  effusion  technique1!2  is  widely  used,  in  varied  guises, 
to  obtain  vapor  and/or  dissociation  pressure  data,  especially  at  high 
temperatures.  In  actual  laboratory  operation  the  conditions  under  which 
effusion  occurs  are  rarely,  if  ever,  the  ideal  conditions  assumed  in 
the  derivation  of  the  simple  Knudsen  equation.  This  report  summarizes 
work  on  several  approaches  designed  to  clarify  the  understanding  of  the 
effusion  process  under  non-ideal  conditions.  More  detailed  introductory 
paragraphs  are  included  in  each  of  the  following  sections. 

SECTION  II 

THEORETICAL  ANALYSIS  OF  MOLECULAR  FLOW 
THROUGH  CONICAL  ORIFICES 

Vapor  pressures  and  the  composition  of  vapors  at  high  temperatures 
are  often  determined  by  effusion  techniques3,  e.g.,  Knudsen  effusion, 
torsion-Kriudsen-effusion-recoil,  and  target-collection  methods;  furthermore, 
mass  spectrometry  of  high  temperature  vapors  is  often  facilitated  by  the 
use  of  an  effusion  cell  as  the  vapor  source3.  Derivation  of  thermodynam¬ 
ically  significant  data  from  the  results  of  effusion  experiments  requires 
knowledge  of  the  relation  between  the  measured  quantities  and  the  equilib¬ 
rium  pressure  in  the  cell;  this  relation  roust  take  into  account,  among 
other  factors,  the  non-ideality4  of  the  geometry  of  physi ca lly-reali zable 
effusion  orifices. 
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The  effect  o£  orifice 


the  rate  of  molccului  effusion 


through  the  orifice  lias  been  analyzed  by  several  investigators5"",  whose 
f-sults  have  usually  been  given  as  numerical  values  for  the  transmission 
probability,  i.e.,  the  probability  that  a  molecule  which  enters  the  orifice 
through  one  end  will  exit  from  the  opposite  end.  The  effect  of  orifice 
geometry  upon  the  angular  distribution  of  molecules  effusing  from  an 
orifice  has  received  less  attention  but  an  analysis10  for  cylindrical 
orifices  has  been  made  and  the  results  have  been  used11  to  obtain  the 
"recoil- force  correction  factors"  applicable  to  torsion- Rnudsen-ef fusion- 
recoil  measurements3. 

In  Part  I12  of  this  report  and  in  a  contribution  to  "Condensation  and 
Evaporation  of  Solids"13  (  subsequently  cited  as  CES)  we  have  presented 
an  analysis  of  the  flow  of  rarefied  gases  through  and  from  conical  orifices. 
Equations  were  derived  for  calculation  of  the  incident  molecular  flu?; 
density  along  the  orifice  walls,  the  transmission  probability  of  the  orifice, 
the  angular  uistribution  of  molecules  leaving  the  orifice,  and  various 
other  functions.  In  Part  I,  but  not  in  CES.  the  results  of  extensive  com¬ 
putations  of  the  various  functions  were  reported;  for  convenience,  one 
table  and  two  Figures  from  Part  1  are  repeated  herein,  and  the  results 
discussed  briefly.  Subsequently,  these  results  are  compared  with  other 
theoretical  results;  comparison  with  experimental  results  is  deferred 
to  Section  III  C  and  IV  of  this  report. 

In  the  last  portion  of  this  Section  we  shall  exploit  the  angular 
distribution  functions  for  conical  oriffees,  which  have  not  previously 
been  available,  to  consider  the  determination  of  an  optimum  orifice  geometry 
(within  t.he  class  of  right  circular  cones)  for  particular  experimental 
configurations. 


A.  RESULTS  OF  THE  NUMERICAL  CALCULATIONS 


In  Figure  1  the  various  critical  orifice  parameters  are  defined  by 
illustration.  Values  of  the  transmission  probability  W,  and  the  recoil- 
force  correction  factor  Jf  are  given  in  Table  1,  and  in  Figures  2  and  J, 
respectively,  for  conical  orifices  with  O^L/jEo^lO  and  with  -90°£Ts4-906, 
i.e,,  for  both  converging  (T  negative)  and  diverging  (T  positive)  conical 
orifices.  Transmission  probabilities  ate  given  only  for  the  diverging 


2 


B 


FLOW 


Figure  1.  Conical  Orifices  and  Parameters:  A,  Diverging 
Configuration,  T  positive;  B,  Converging  Configuration, 
T  negative. 


5 


TABLE  I 

TRANSMISSION  PROBABILITIES  AND  RECOIL  FORCE  CORRECTIONS 
FOR  CONICAL  ORIFICES 


(L/r  )fl 

W 

i 

(fVrm2)a 

0,1 

0,952399 

0.968322 

0.2 

0.909215 

0.937308 

0.4 

0.82408 

0.87847 

0.6 

0.77115 

0.82471 

0,8 

0.71778 

0.77620 

1 .0 

0,67198 

0.73269 

2,0 

0.5142 

0.5725 

4.0 

0.3566 

0.4024 

6.0 

0.2754 

0.3125 

8.0 

0.2253 

0.2564 

10.0 

0.1909 

0,2177 

0.1 

0.967347 

0.986104 

0,980330 

0.2 

0,93835 

0.97200 

0.96127 

0,4 

0.88938 

0.94451 

0.92584 

0.6 

0.84990 

0.91904 

0.89445 

0.8 

0.8176 

0.8960 

0.8670 

1.0 

0.7908 

0.8756 

0.8432 

2,0 

0.7058 

0.8035 

0,7621 

4.0 

0.6348 

0 , 7381 

0,6899 

6.0 

0.6051 

0.7107 

0.6589 

8.0 

0.5895 

0.6970 

0.6425 

10.0 

0.5802 

0.6892 

0.6328 

0.1 

0.978646 

0.999546 

0.988540 

0.2 

0.96027 

0.99816 

0.97764 

0.4 

0.93057 

0.99407 

0.95804 

0.6 

0.90793 

0,98941 

0.94150 

0.8 

0.89034 

0.98492 

0.92770 

1.0 

0.87642 

0.98087 

0.91625 

2.0 

0.8370 

0.9682 

0.8813 

4.0 

0.8108 

0.9633 

0.8564 

6.0 

0.8022 

0.9660 

0.8480 

8.0 

0.7984 

0.9699 

0.8442 

10,0 

0.7963 

0.9736 

0.8421 
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<f/r  *) 


0.98591b 

0.976151 

0.95973 

0.94812 

0.9396 

0,93338 

0.9177 

0.9095 

0.507  3 

0.9065 

0.9060 

0.992680 

0.987008 

0.97902 

0.97389 

0.97044 

0.96806 

0.96288 

0.9607 

0.9603 

0.9601 

0.9599 

0.994775 

0.990881 

0.985662 

0.98249 

0.98047 

t\  *\  ->  i  «-» 

KJ  ,  'J  /  J  l  / 

0.9764 

0.9754 

0.9752 

0.9751 

0.9750 

0.996420 

0.993877 

0.99066 

0.98883 

0.98772 

0.98701 

0.98570 

0.9853 

0.9832 

0.9641 

0.98  -1 


1.008954 
1 .016307 
1 .02782 
1  .03652 
1 .0434 
1  .04919 
1.0687 
1.0908 
1  .1044 
1.1140 
1  .1212 

1  .014589 
1  .026941 
1  .04683 
1  .06220 
1  .07444 
l  .08446 
1  .11632 
1.1462 
1.1618 
1.1720 
1  . 1  793 

1 .016077 
1 .029628 
1 .051287 
1.06782 
1 .08086 
1 .09143 
1  .1242 
1  .1 530 
1  .1684 
1.1778 
1 .1844 

1 .016729 
1 . 030684 
1 .05267 
1  .06918 
1 .08202 
1 .04231 
1 .12355 
1 . 1 504 
1  .  1  636 
1.1717 
1.17  74 


0.993868 

0.988198 

0.97848 

0.97079 

0.9648 

0.96004 

0.9473 

0.9401 

0.9381 

0.9373 

0.9369 

0.997083 

0.994494 

0.99034 

0,98732 

0.9851  3 

0.98352 

0.97974 

0.9780 

0.9776 

0.9775 

0,9  7  74 

0.998104 

0.996468 

0.993949 

0.99221 

0.99101 

il.vvOib 

0.9883 

0.9875 

0.9874 

0.9873 

0 . 987  3 

0.99883! 

0.497855 

0.99642 

0.49549 

0 . 44489 

0.99447 

0.99363 

0.9933 

0 , 4933 

0.4932 

0.9932 


T 


(L/r  ) 1 
m 

1  j 

£ 

i£*/r 

ni 

0,1 

0.988591 

1.015712 

0.999644 

0.2 

0.997715 

I .028324 

0.999374 

0  .4 

0.996761 

1 .0^7307 

0.999026 

0.6 

0.996304 

1.060889 

0.998831 

0.8 

0. 9960o 

1.07108 

0.99673 

1 .0 

0.99592 

1 .07904 

0.69866 

2.0 

0.99571 

1 . 10207 

0.99854 

4.0 

0.9957 

1 . 1205 

0.9985 

6.0 

0.9957 

1 . 1289 

0.9985 

8.0 

0.9956 

1.1338 

0.9985 

10,0 

0.9956 

1 . 1372 

0.9985 

0.1 

0.999631 

1.011 961 

0.999936 

0.2 

0.999450 

1 .020864 

0.999895 

0.4 

0.999296 

1.033229 

0.999853 

0.6 

0.99924 

1.04141 

0.99983 

0.8 

0.99921 

1,04724 

0.99982 

1  .0 

0.99920 

1 .05160 

0,99982 

2.0 

0.9992 

1.0635 

0.9998 

4.0 

0.9992 

1.0721 

0.9998 

0.1 

0. 999960 

1 .006048 

0.999999 

0.2 

0.999957 

1 .009755 

0.999997 

0.4 

0.999951 

1 .014072 

0.999996 

0.6 

0.999  9  5 

1 .01651 

1.00000 

0.8 

0.99995 

1 .01809 

1  .00000 

1.0 

0.99995 

1 .01919 

1 .00000 

2.0 

1  .004)0 

1.0219 

1 .0000 

4.0 

1.0000 

1.0236 

1  .0000 

r  «  radius  of  smaller  end  of  orifice.  g)  is  the  appropriate 

“in 

recoil  force  correction  for  Converelne  orifices  when  (h/r)  for  the 
orifice  is  calculated  with  r  rather  than  rG  (see  Figure  1  and 

—iu  - 

references  12  and  lit)* 


Log  (L/rc) 


Figure  2.  The  Transmission  1’iobability  W  as  e  fu 
log  (h/r0)  at  various  values  of  T. 


configuration  because,  for  a  given  orifice  geometry,  U  is  the  same  fox  gas 
flow  In  either  direct  ion*'* 14.  The  logarithmic  scale  on  the  abscissa  of 
Figures  2  and  3  is  used  only  to  simplify  presentation  of  the  results. 

Table  1  and  Figure  2  reveal  no  unexpected  variation  of  W;  it  decreases 
with  increasing  L/'ro  and  increases  as  |T|  increases,  as  would  be  predicted 
intuitively.  However,  Figure  3  reveals  an  unexpected  variation  of  _f  wfth 
T:  for  _T  >  22°,  t  is  great ex  than  unit v.  Hcnrc,  the  recoil  force  generated 

by  molecular  effusion  from  a  conical  orifice  with  T  >  22°  is  greater  than 
that  from  the  corresponding4  Ideal  orifice.  This  result  arises  because 
(1)  the  overage  angle  at  which  molecules  efiuse  from  conical  orifices  is 
smaller  than  for  the  ideal  orifice,  and  (,°)  as  T  Increases  above  20°,  the 
transmission  probability  rapidly  approaches  unity;  the  moniental  component 
along  the  conical  orifice  axis  is  therefore  greater  than  that  along  the 
ideal  orifice  axis  and  the  recoil  force  fs  greater.  This  ‘’focusing" 
eifecL  of  conical  orifices  is  a  maximum  at  the  precise  value  depending 

on  (b/xo). 

Numerical  values  for  various  other  quantities,  c.g,,  the  incident 
molecular  flux  on  the  orifice  wall  and  t he  angular  distribution  of  effusing 
molecules,  have  been  calculated  but  arc  not  presented  here;  an  extensive 
tabulation  o£  these  functions  is  available  elsewhere 14 . 

B.  COMPARISON  WITH  Or  HER  THEORETICAL  RESULTS 

It  was  shown  in  Part  1  and  in  CES  that,  where  comparison  was  possible, 
the  results  presented  were  in  agreement  with  those  of  other  workers.  In 
particular,  our  results  are  in  excellent  agreement  with  the  transmission 
probabilities  ol  lezkowski,  et  a]."  for  conical  orifices,  and  with  the 
closely-bounded  values  for  the  transmission  probabilities  of  cylindrical 
orifices  obtained  by  DeMarcus  and  Hopperf’.  In  view  of  the  different 
approaches,  the  independent  verification  ol  derivations,  and  the  Tonsider- 
ably  different  numerical  techniques  used  i  i  solving  the  integral  equations, 
by  four  groups5; 14  there  can  be  little  doubt  of  the  valfdity  of 
the  numerical  values  for  the  transmit slon  probability  for  cylindrical  and 
conical  orifices;  these  values  apply  rigorously  only  to  the  assumed  model, 
of  course. 


We  know  of  no  other  analysis  of  annular  distribution  lot  conical 
orifices,  and  for  cylindrical  orifices  only  of  Clausing^  analysis10. 
Clausing  assumed  that  the  normalized  incident  density  (or  molecular  flux) 
on  the  walls  of  a  cylindrical  orifice  could  be  expressed  as  a  linear 
function  of  the  distance,  along  the  orifice  axis,  from  the  entrance. 

The  accuracy  of  this  approximation  has  been  discussed  by  Edwards16  who 
has  also  shown  that,  with  this  approximation,  Clausing^  distribution 
function  i°,  tl  t(0) 

may  be  derived  in  detail  directly  irom  our  more  general 
distribution  tunctions  which  apply  to  conical  and  cylindrical  orifices; 
it  has  also  been  possible  to  integrate  in  closed  form  the  resulting 
angular  distribution  equation  for  cylindrical  orifices. 

One  further  argument  may  be  advanced  for  the  validity  of  our  angular 
distribution  results.  In  solving  the  various  equations 13  the  first  quantity 
obtained  is  i'(x)  the  normalized  incident  density  on  the  orifice  walls; 
these  values  of  V(x)  are  then  used  to  obtain  W  in  a  direct  way.  The  same 
values  of  i|<(x)  are  used  in  the  more  involved  computation  of  W  by  numerical 
integration  of  the  angular  distribution  functions.  The  two  values  of  W 
are  in  good  agreement  with  each  other13  and,  ns  already  noted,  with  the 
values  of  Iczkowski,  et  al.8,  which  would  seen  to  indicate  that  both  i|i(x) 
and  the  angular  distribution  functions  are  correct. 

It  is  difficult  to  compare  our  results  with  those  of  Davis,  at 
their  paper  lacks  detail  in  both  derivation  and  results.  Our  impression 
is  that  their  derivation  is  not  rigorous.  In  any  case  the  result.?  shown 
in  their  Figure  7  for  the  transmission  probability  of  ’’conical  nozzles" 
appear  to  be  higher  than  those  of  Table  1  by  several  percent  for  smell 
angles  T  and  lower  for  large  angles  T. 

Sparrow  and  Jons son"3  have  also  analyzed  the  mass  flow  through  conical 
orifices  and  have  used  the  results  to  formulate  analysis  of  energy  transfer 
between  gas  end  orifice  walls.  If,  in  their  equation  ( 12)  and  the  ordinate 
of  their  Figure  3,  Lz  fs  made  zero,  the  resulting  cerro  is  equivalent  to 
our  transmission  probability  W.  The  results  in  their  Figure  p  appear  to 
be.  in  good  agreement  with  values  of  W  in  Table  is.  and  f ;(x/T^  of  their 
figure  (2)  appear  to  agree  with  out  Sparrow  and  Joncaon  obtain 

a  number  of  useful  but  seldom  used  relations  ( thelj  equations  some 

of  which  we  have  also  generated13.*  l4. 
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Richley  and  colleagues  liav.  recently  reported  a  sctics  of  analyses 
of  molecular  flow  through  cylindrical  tubes13^ c,  converging  and  diverging 
tubes3"1-1  and  slotsl!;>b>d>  and  through  cylindrical  tubes  '..•ill;  contribution 
from  surface  diffusion iyc»M  .  The  investigation  by  Cook  and  Rich Icy  of 
angular  distribution  from  cylindrical  orifices  is  based  on  Clausing  s 
analysis10  which  has  already  been  correlated  with  the  present  results. 

In  their  analysis19^  of  flow  through  converging  and  diverging  tubes  and 
slots,  Richley  and  Reynolds  obtain,  by  iterative  solution  of  the  appro¬ 
priate  Fredholm  integral  equation,  values  for  the  normalized  incident 
density  y(x)  on  the  orifice  wall  (their  "flux  ratio  ]" )  and  use  these 

values  to  obtain  the  flux  distribution  over  the  plane  of  I. be  exit  end  of 
the  orifice;  finally,  the  exit  plane  flux  distribution  is  integrated  over 
ttie  exit  area  of  the  orifice  to  obtain  the  transmission  probability  W 
(their  Ft)’  The  minor  discrepancies  which  exist  between  their  results 
and  ou»>;  appear  to  arise  from  ( l)  their  use  of  an  iterative  solution  for 
1  rather  than  cut  more  direct  solution13?  14  for  i}- ( x ) ,  and  (2)  their 
introduction  of  an  additional  numerical  integration  (i.e.,  to  evaluate 
the  exit  plane  flux  distribution)  in  the  computational  sequence  leading 
to  W,  rather  than  calculation  of  W  directly  from  the  i['(x)  In  one  step13. 


C.  OPTIMIZATION  OF  ORIFICE  GEOMETRY 


An  orifice  geometry  will  be  considered  optimum  if,  for  any  given  race 
of  effusion,  the  quantity  being  measured  is  maximized  with  respect  to  the 
orifice  parameters  T  and  ( L/'ji_q) ;  the  effusion  rate  may  always  be  adjusted, 
if  necessary,  by  varying  the  orifice  area  a 0  while  maintaining  a  fixed 
geometry,  i.e.,  fixed  values  of  T  and  ( b/ro)  •  This  criterion  for  optimi¬ 
zation  is  directly  applicable  when  any  one  of  the  following  typical  con¬ 
ditions  exist:  (a)  a  very  small  amount  of  sample  is  available;  (b)  intro¬ 
duction  of  new  sample  into  the  apparatus  requires  a  comparatively  long 
down-time  (as  in  mass  spectrometry);  (c.)  the  material  under  study  vaporizes 
incongruent ly  and  vaporization  characteristics  change  with  the  composition 
of  the  condensed  phase;  and  ( d)  a  low  evaporation  coefficient  imposes  the 
need  to  minimize  the  total  flux  from  the  cell,  thereby  minimizing  the  dis- 
- la cement  of  the  actual  pressure  from  the  equilibrium  value.  Under  each  of 


these  conditions  it  is  clearly  desirable  to  maximize  the  measured  quantity 
while  simultaneously  minimizing  the  total  rate  of  effusion  of  sample  from 
the  cell. 

We  shall  now  determine  the*  orifice  geometry  which  maximizes  each  of 
the  following:  (  l)  the  near-axial  flux  density  of  the  effusing  molecular 
stream;  (2)  the  recoil  force  generated  on  the  effusion  cell;  and  (3)  the 
force  exerted  on  a  target  suspended  in  the  effusing  gas  stream. 


1.  Optimization  o.f  Molecular  Beam  Intensity  on  end  near  the  Orifice 
Axi s .  In  Tart  I12  and  CES*3  we  obtained  expressions  for  what  we  shall 
here  call  t lie  incremental  transmission  probability  _dWg,  ..e.,  t.he 
probability  that  a  molecule  which  enters  one  end  of  the  orifice  will 
exit  from  the  opposite  end  into  the  incremental  solid  angle  dm  located 
at  angle  0  from  the  orifice  axis.  With  the  solid  angle  Ac  expressed  in 
spherical  coordinates  and  the  assumption  of  circular  symmetry  in  the 
distribution  around  the  orifice  axis,  the  equations  may  be  written  as 

dWe  =  2Qn  SinQ  Cos0  d9;  ( l) 


the  _Qn  ( n  =  1,2,3)  are  complicated  functions  of  the  orifice  parameters 
T,  (jb/xo) >  and,  depending  on  the.  range  in  which  0  lies,  also  of  0. 

The  integral  ol'  dWn  over  0  i  3  i  tt/2  is  just  the  transmission  probability 
W  ( designated  Wg  in  Part  I  and  CES  to  distinguish  a  calculation  from 
angular  distribution  considerations). 

We  now  define  th;  quantity  Wy, 


v  *y 

Wy  ■»  Pq  dWg  =>  Pq  2Qn  Sing  Cose  d0. 


(2) 


which  may  be  called  the  angular  transmission  probability,  i.e.,  it  is 
the  ratio  of  the  number  Ny  of  molecules  which  effuse  per  second  into 
a  cone  which  is  coaxial  with  the  orifice  and  has  semi-apex  angle  Y, 
to  the  number  N0  of  molecules  which  enter  the  orifice  through _rQ  per 


second. 


The  quantity  to  be  maximized  is  the  ratio  of  NY  to  the  total 
number  Ng  of  molecules  which  effuse  from  the  orifice,  i.e.,  Ny/N,  ** 
(_Ny/JSf0) °)  "  ~y/—>  and  irow  equation  (2), 


WY/W  -  (2/W)  jJJ  Qn 
If  we  restrict  Y  to  Range  I13'  14 


SinG  CosG  dQ. 
(i.e,,  to  0  s  V  i 


OJ 


t]),  5n  becomes  <Ja 
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and  is  independent  of  0;  equation  (3)  may  be  now  integrated  to  obtain 
Wy/W  =  (Qa/W)  Sin"  Y;  0  s  V  i  |T|,  (4) 

Quite  obviously,  the  fraction  of  effusing  molecules  which  flow  into 
the  cone  defined  by  the  orifice  axis  and  the  angle  y  can  be  varied  by 
varying  y;  the  pertinent  problem  is  to  maximize  that  fraction  for  a 
ven  y  determined  by  a  particular  experimental  apparatus.  Considering 
then  that  Y  is  fixed,  we  re-write  equation  (4)  as 

Wy/(W  Sin2  Y)  =  (Q^W),  (OsKslIl),  (5) 

and  note  that  the  quantity  (.Qi/w)  now  to  be  maximized  is  a  function 
solely  of  orifice  geometry. 

A  plot,  prepared  from  data  tabulated  elsewhere14,  of  (jji/w)  vs.  T 
for  various  (h/r0)  is  given  in  Figure  4.  It  is  readily  apparent  that 
a  maximum  exists  at  T  =  0°  for  each  (l/to)  and  that  the 
maximum  becomes  sharper  as  (L/ro)  increases.  He  now  note  that  under 
the  restriction  0  £  Y  s;  |T|,  when  T  =  0°,  Y  can  have  only  the  value 
0° ;  hence,  the  points  in  Figure  4  for  T  =  0°  apply  only  to  the  beam 
intensity  alone,  the  orifice  axis  where  Y  =  0°. 

To  determine  whether  Wy/(W  Sin'  Y)  is  also  a  maximum  at  T  =  0° 
when  Y  has  typical  experimental  values  of  5°  and  103,  values  of  Wy 
for  cylindrical  orifices  ( T  =  0° )  have  been  calculated,  in  effect  by 
integration  of  equation  (2),  but  in  fact  from  a  closed” form  solution14* 15. 
The  results,  plotted  as  WY/(W  Sin”’y),  are  in  Figure  5  in  which  the 
curves  for  T  s  -10”  and  for  T  2  +10°  duplicate  those  of  Figure  4 
since  in  these  ranges  (y  i  |T|). 

It  Is  evident  in  Figure  5  that  for  y  =  5%  W-./(W  Sin2  Y)  is 
larger  for  T  «*  0°  than  for  any  T  2  r/  ;  similarly,  for  y  *»  10° 

(WY/W  Sin2  y)  is  larger  for  T  =  0°  than  for  any  T  >  10°.  It  appears 
to  be  true  in  general  (although,  we.  have  not  made  extensive  computations 

*  - - -,r  I -  »•  I,  r  ~  (  T  /  _  \ - 1  _  „  J - /  TT  /  T.7  OJ-2  *w\ 

tv  twiitinu  XrHJL.9j  tnac  wi  a  \±i/ JiO/  cauu  a  ^ivcu  *  }  VJ£y  f  XL  ,JJ-“  T/ 

is  larger  for  T  =  0°  than  is  (,2i/W)  for  any  T  s  y. 

For  a  short  orifice  (L/ro  <  2.0^  and  a  given  y,  the  two  quantities 
(SV^)x=  anc*  (W-y/bT  Sin2  y)q=o  are  trie  same  within  1$;  for  longer 
orifices  the  difference  is  somewhat  greater. 
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There  remains  the  (jucouiua  u I  the  relative 
for  T  «  0°  and  for  some  T  <  y,  c.g.,  for  T  13  0° 
For  these  cases  (W  /W  SlnH  y)  is  larger  tor  the 
(T  >  0)  and  Increase*  slowly  with  increasing  T. 


value  of  (Wy/W  Sin'"  y) 
and  T  «  5°  with  Y  **  10°. 
conical  orifice 


2.  Optimization  of  the  E f fusive  Rccoi 1  Force 

In  CES 13  (equation  19 )  it  was  demonstrated  that  the  "incremental 
recoil-force  correction  factor"  _dfg  ( designated  dFy  in  CES).  (l.e., 
the  recoil  force  dfl  genera 1 1  1  on  a  cell  by  effusion  of  molecules  into 
the  incremental  solid  angle  _dju  at  angle  6  from  the  orifice  axis, 
divided  by  the  total  recoil  force  0^  for  the  corresponding4  ideal 
orifice:  dig  =  is  related  to  _dWg  by 


df g  =  (3/2)  Cose  dW0 
or,  with  equation  (l), 

dfg  <=  3Qn  SinS  Cos2  G  d0.  (6) 

The  (tatad;)  recoil  force  <cx«rreotion  factor  11> 3  is  then  given  :by 
e=rt/  2  3  bn 

f  13  Je=0  dfQ  a  nil  Jan  5Q"  Sin0  Cos2  0  d0;  (7) 

The  integration  limits  ( anj  bn)  are  discussed  in  CES  and  reference  14. 

The  factor  _f  is,  of  course,  the  quantity  tabulated  In  Table  1. 

Optimum  orifice  geometry  in  tills  situation  requires  maximization 
of  recoil  force  for  any  given  rate  of  effusion,  which  is  equivalent 
to  maximizing  the  ratio  (_f/W).  In  Figure  6  this  ratio  is  plotted  vs. 

T  fox  various  (h/i.o)-  xl  is  eviaenc  from  Figures  2  and  5  that  a  maximum 
must  occur  in  a  plot  of  ( jE/W)  vs.  T.  However,  it  is  rather  surprising 
that  this  maximum  occurs  at,  or  very  near,  T  =  +30°  for  a  very  wide 
range  of  h/r0.  It  is  quite  clear  from  Figure  6  that  the  optimum 
orifice  for  recoil  force  measurements  would  have  a  semi-apex  angle  T 
of  +30°  (i.e.,  diverging)  and  would  be  as  long  as  practicable  (within 
the  range  of  jJx.o  covered  in  Figure  6). 


3*  Optimization  of  the  Force  Exerted  on  a  Target  in  _an  Effusing 
ftolecular  Beam:  We  shall  restrict  our  consideration  to  cases  In 
which  the  target  la  circular  and  coaxial  with  the  orifice,  and  In 
which  the  molecules  striking  the  target  either  all  condense  or  all 
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revalorize.  The  angle  subtended  at  the  orifice  by  the  target  is 
designated  21;  V  then  has  the  same  meaning  as  in  equation  2. 

If  all  molecules  condense  on  the  target,  the  ratio  of  the  force 
Fy  exerted  on  the  target  to  the  force  F0  which  would  be  exerted  by 
all  molecules  effusing  from  the  corresponding4  ideal  orifice,  is  given 
by  equation  (7)  with  integration  over  0  i  6  s  y; 

fy  -  Fy/F0  -  Jg  3Qn  SinG  Cos2e  do.  (8) 

If  Y  is  restricted  to  0  <;  V  £  |T  j  so  that  only  JJi  (*•  constant)  is 
required  in  the  integration,  equation  (8)  becomes 

fY  ■=*  Qi(  1  -  Cos3y) . 

We  wish  to  maximize  this  force  at  a  given  flow  rate  and  for  a  given  Y; 
therefore,  we  write 


The  quantity  to  be  maximized  13  (Jj/W,  as  in  subsection  1  above,  and 
the  arguments  given  there  apply. 

If  the  molecules  revaporize  from  the  target,  there  will  be 
exerted  on  the  target  an  additional  force  which,  with  all  extra¬ 
orifice  parameters  fixed,  will  depend  on  the  molecular  flux  onto  the 
target3.  Maximization  of  this  flux  at  a  given  total  flux  from  the 
orifice  was  the  subject  of  subsection  1  above;  again,  Qi/W  is  the 
quantify  to  be  maximized. 


SECTION  III 


MEASUREMENT  OF  ANGULAR 
DISTRIBUTION  OF  MOLECULAR  FLOW 
THROUGH  CONICAL  ORIFICES 

Of  the  various  quantities  which  are  derived  in  the  theoretical  anal¬ 
ysis  12 1  14  (see  Section  II)  and  which  are  amenable  to  experimental  study, 

the  most  critical  is  the  angular  distribution  of  molecules  effusing  from 
an  orifice,  i.e.,  the  variation  with  0  of  the  number  dN  of  molecules  which 
pass  per  second  from  the  orifice  into  the  incremental  solid  angle  da)  located 
at  angle  6  from  the  orifice  axis.  Tire  theoretical  analysis  predicts  that 
dN^  is  proportional  to  cos  0;  J2n  is  <-l'e  complicated  function  of  orifice 
parameters  which  arises  from  the  non-ideality  of  the  orifice  (for  the  ideal 
orifice,  (Jn  is  always  unity).  Measurement  of  dN  or  an  equivalent  quantity 
would  provide  experimental  data  which  could  be  compared  directly  with  the¬ 
oretical  values  for  £n  cos  0. 

A.  EXPERIMENTAL  APPARATUS 

To  accomplish  these  measurements  the  apparatus  described  briefly  in 
this  section  (and  in  detail  in  Part  II*21)  has  been  constructed.  It  is 
designed  to  allow  a  study  of  the  ef fusion  of  a  permanent  gas  (e.g.,  Na, 

He,  COa)  at  any  suitable  pressure  from  any  orifice  with  a  geometry  which 
can  be  machined  into  a  small  circular  plate.  Permanent  gases  are  used 
as  effusants  so  that  the  apparatus  can  he  operated  at  room  temperature; 
concern  that  the  reservoir-orifice  system  may  not  be  isothermal  is  there¬ 
by  minimized.  For  this  advantage  the  ability  to  study  the  effusion  process 
as  a  function  of  temperature  is  sacrificed. 

The  angular  distribution  of  effusing  molecules  is  determined  by  a 
molecular  beam  method22  incorporating  a  modulated  beam  technique23.  The 
reservoir  from  which  the  molecules  effuse  can  be  rotated  on  an  axis  which 
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passes  through,  and  ts  parallel  to,  the  outer  face  of  the  orifice  under 
study.  Two  stationary  collimating  orifices  and  the  effusion  orifice  under 
study  define  a  molecular  beam,  the  beaiu  is  modulated  by  a  mechanical  chopper, 
and  its  intensity  is  determined  by  a  neutral-beam  detector. 

The  apparatus,  a  diagrammatic  horizontal  cross-section  of  which  is 
shown  in  Figure  7,  consists  of  five  principal  components:  ( l)  The  Main 
Vacuum  Chamber,  which  can  be  maintained  at  a  pressure  very  low  with  respect 
to  the  pressure  in  the  gas  reservoir,  and  wnich  contains  the  rotating 
eifusing  cell;  ( 2)  The  Gas  Kcscrvolr.  a  large  chamber  from  which,  gas 
flows  to  the  effusion  cell  and  in  which  the  pressure  can  be  kept  constant; 

(3)  The  buffer  Chamber,  a  small  independently-pumped  chamber  which  is 
separated  from  the  main  chamber  by  a  plate  containing  the  first  collimating 
orifice,  and  which  contains  a  chopper  capable  of  interrupting  the  beam 
about  one  hundred  times  per  second;  (h)  The  Detector  Chamber,  an  independ¬ 
ently-pumped  volume  which  is  separated  from  the  buffer  chamber  by  a  plate 
containing  the  second  collimating  orifice,  and  which  contains  an  electron- 
impact  molecular  beam  detector;  (5)  The  Detector  Electronics,  which  con¬ 
sists  of  a  power  supply  for  the  beam  ionizer  and  a  system  to  amplify, 
measure,  and  record  the  ion  current  from  the  ion  collector. 

Each  of  these  components  was  described  in  detail  in  Part  IIal;  we  shall 
describe  here  only  significant  modifications. 

1.  Modulation  of  Mocular  Ream:  As  experimental  techniques  were 
refined,  it  became  apparent  that  the  frequency  of  the  mechanical 
beam  chopper  (Figure  ‘p,  Part  II)  was  not  stable;  the  cause  was 
found  to  be  an  overloaded  driving  motor.  When  the  motor  was  re¬ 
placed  by  one  with  higher  torque,  a  new  motor  mounting  and  chopper 
housing  (items  6  ane  12,  Figure  5,  Part  II)  was  fabricated  from 
stainless  steel  (rathei  than  brass).  The  new  drive  mechanism  per¬ 
formed  satisfactorily  at  times,  but  at  others  exhibited  erratic 
chopping  frequency.  This  difficulty  arose  from  binding  in  the 
bearings  which  support  the  chopper  shaft;  the  binding  in  turn 
was  caused  by  a  slight  misalignment  of  the  new  motor  mount/chopper 
housing  when  it  was  welded  into  the  vacuum  system. 

At  this  point  an  electronically-driven  tuning  fork  chopper 
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The  Main  Vacuum  Chamber 
The  Cas  Reservoir 
The  Fuller  Chamber 
The  Detector  Chamber 
The  Rotating  El  fusion  Cell 
The  Ream  Collimating  Orifices 
The  Beam  Chopper 
The  Beam  Ionizer 

Connectors  lor  Tygon  lubulation  Which  Carries  Gas 
from  t  lie  Reservoir  to  the  Rotating  Cell 
Tumping  "l"s  Which  Accommodate  hi  quid  Nitrogen  Traps 
Glass  Windows 

Rotary  Vacuum  Seal  Through  Which  the  Cell  is  Rotated 
Brass  Bellows 
Globe  Valve 

Copper  Pipe,  1  1/-  -  in.  Diameter 
Port  for  Attaching  ^Ij-Liter  Stainless  Steel  Tank 
Valve  for  Introducing  Eifusant  Gas 
Connections  to  the  Lquibar  Pressure  Metei 


Figure  7  -  Diagram  of  the  experimental  Apparatus. 
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(type  I4O)  was  obtained.  Nominal  frequency  of  the  chopper  is  170  Hz; 
the  chopping  vanes  attached  to  the  tines  are  10  mm  high  and  have  a 
maximum  aperture  ot  8-10  mm.  An  electronic  signal  synchronized  with 
the  mechanical  oscillations  of  the  tuning  fork  is  an  inherent  feature 
of  the  driving  circuit.,  and  is  readily  available  for  use  as  reference 
signal  to  the  lock-in  amplifier. 

The  new  chopper  was  mounted  _in  the  BUFFER  CHAMBER  (Figure  7) , 
but  directly  on  the  flange  which  is  welded  onto  the  MAIN  VACUUM  CHAMBER 
and  which  mates  with  the  BUFFER  CHAMBER  flange.  A  new  BUFFER  CHAMBER 
equipped  with  aluminum- foi 1-sealcd  flanges'  was  fabricated  fiom 
stainless  steel. 

No  particular  problems  have  been  encountered  with  the  toning 
fork  chopper;  it  was  used  in  obtaining  essentially  ;  the  data 

reported  in  Section  B. 


o.  Detector  Chamber :  The  Detector  Chamber  described  in  Part  11 
operated  satisfactorily,  except  that  the  ultimate  pressure  was  rel¬ 
atively  high  (~10"c  torr).  Consequently,  the  life  of  the  cathode 
was  seriously  shortened  and  the  emission  current  available  from  the 
cathode  was  low.  A  new  chamber  was  fabricated  from  stainless  steel; 
aluminum  foil  flanges-'M  were  used,  and  all  elastomer  0* rings  were 
eliminated;  a  more  efficient  liquid  nitrogen  trap  ( "Cryosc-rb")  was 
inserted  between  the  2"  diffusion  pump  and  the  chamber. 
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reaches  >xlO_B  torr  and  operates  at  pxlO’7  tore  with  a  beam,  from 
a  0.1  torr  source,  entering  the  chamber. 

The  basic  beam  detector  design  and  electronics  have  undergone 
no  significant  change.  In  the  improved  vacuum  system  the  beam 
defector  has  performed  quite  satisfactorily;  cathode  liie  is  remark¬ 
ably  lengthened  (no  quantitative  data,  as  yet)  and  the  emission 
current  is  increased  by  a  factor  of  live. 


3.  Valve?  Tire  sensitivity  of  the  beam  detector  cathode  to  oxygen 
mode  it  desirable  to  keep  the  detector  chamber  evacuated  continuously. 
To  achieve  isolation  of  the  Detector  (and  Bui  for/  Chamber  virile  the 
Main  Chamber  was  opened,  c. g. ,  to  change  orifices,  a  sliding. 
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sealed  valve  was  installed  in  the  Main  Chamber  in  such  a  way  that 
closure  of  the  valve  sealed  the  first  collimating  orifice  from  Lite 


Main  Chamber.  The  design  of  the  valve  was  adapted  from  that  of 
Sheffield1^.  It  has  worked  satisfactorily. 

ll .  Alignment  of  Orifices,  The  apparatus  was  designed  and  constructed 
to  permit  direct,  visual  alignment  of  the  beam  defining  orifices, 
if  the  Faraday  ion  collector  is  removed  from  the  beam  detector 
assembly.  Alignment  of  the  effusion  orifice  and  the  various  colli¬ 
mating  orifices  is  rather  easily  accomplished  by  viewing  with  a  low- 
powered  telescope  along  the  beam  axis.  However,  with  only  this  check 
on  alignment  it  is  possible  for  the  plane  of  a  given  orifice  to  be 
tilted  appreciably  from  the  desired  perpendicularity  to  the  beam 
axis  and  the  till,  be  undetectable  through  the  telescope. 

Therefore,  to  insure  parallelism  of  the  planes  of  all  beam- 
defining  orifices  while  simultaneously  establishing  all  orifices  to 
be  coaxial,  a  second  alignment  procedure  was  adopted.  A  low-powered 
telescope  was  fitted  with  a  Gaussian  eyepiece  .and  the  telescope 
operated  as  an  autocollimator;  A  light  source  on  the  side  of  the 
eyepiece  illuminated  a  cross-hair  the  image  of  which  was  then  pro¬ 
jected  through  the  telescope  onto  an  optically  fl.st  (both  sides) 
front-surfaced  mirror  attached  firmly  to  the  orifice  plate  being 
aligned;  the  orifice  plate  was  adjusted  until  the  reflected  image 
of  the  cross-hair  coincided  in  the  eyepiece  with  the  image  of  the 
actual  tuoss-iiair;  coincidence  oi  the  two  images  requires  that  the 
mirror  (orifice  plate)  be  perpendicular  to  the  light  (molecular) 
beam. 


B.  EXPERIMENTAL  RESULTS 


Given  the  apparatus  described  above  and  the  desired  experimental 
measurement,  i.e,,  the  angular  distribution  of  the  effusing  beam,  experimental 
procedure  was  rather  straightforward.  After  a  given  orifice  was  inserted 
into  the  Rotating  Effusion  Cell  and  aligned,  the  entire  system  was  evacuated, 
e.g.,  overnight.  At  the  beginning  of  a  run  the  Gas  Reservoir  was  isolated 
from  the  vacuum  pumps  (by  closing  valve  14 .  Figure  7)  -nd  filled  with  gas 
(helium,  nitrogen,  nitrous  oxide)  to  the  desired  pressure  as  measured  by 
a  variable  capacitance  sensor  (Equibar  120).  Cases  were  obtained  from  high- 
pressure  cylinders;  pressure  in  the  reservoir  was  controlled  with  a  variable- 
leak  valve  (Type  9101-M)  in  the  line  between  the  usual  cylinder  regulator  and 
the  reservoir.  After  steady  state  flow  was  established  throughout  the  gas 
flow  system,  the  pressure  in  the  reservoir  remained  surprisingly  constant 
(i.1  to  2%);  a  precision  pressure  regulator  originally  planned  for  insertion 
in  the  gas  line  between  the  cylinder  regulator  and  the  variable- leak  valve 
was  not  required.  I’urity  of  the  helium  and  nitrogen  used  was  >99.5 %  and 
of  the  N20  >98.0$. 

While  steady  gas  flow  was  being  established,  the  electronic  circuitry 
was  energized,  adjustment  of  the  lock-in  amplifier  checked,  and,  in  partic¬ 
ular,  the  filament  current  of  the  beam  detector  was  adjusted  to  provide  an 
electron  emission  current  of  10mA,  The  isolating  valve  (section  III. A. 3) 
was  then  opened;  the  molecular  beam  could  then  pass  to  the  detector  and 


measurements  were  begun. 

With  the  gas  pressure  in  the  reservoir  constant,  beam  intensities  in 
arbitrary  units  (i.e.,  the  output  from  the  lock-in  amplifier)  were  recorded 
with  the  Effusion  Cell  rotated  to  orientations  varying  by  angular  increments 
of  5°  between  0°  and  +90°  and  also  between  0°  and  -90° .  Typical  concordance 
between  data  for  +@  and  for  -0  is  illustrated  in  Figure  8. 

Angular  distribution  data  have  been  obtained  for  five  orifices  with 


three  gases  et  various  pressures; 


parameters. 

In  Figures  9~l4  the  results  of  various  measurements  are  plotted  as  Ig 
vs.  0;  _Ig  is  the  relative  molecular  beam  intensity  normalized  to  1,0  at 
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TABLE  2 


ORIFICES,  GASES,  AND  PRESSURES 
USED  IN 

ANGULAR  DISTRIBUTION  STUDIES 

- -  - - EEMSl-lRE^^^jj^torr 


Orifice 


Number 

T,deg 

L/rm 

Lm*  tni 

GAS 

lie 

Na 

N20 

1 

23.05 

4.010 

0. 0460 

190.,  70.0, 
45.0,  20.2 

070.,  200., 

72.7,  4i.o, 
27.5,  20.2, 

7.0 

195.,  69.0, 
45-0,  22.0, 
6.9 

2 

8.50 

10.08 

0.0455 

190.,  71.0, 

41.5,  21.0 

660.,  210., 
70.0,  22.0, 
7.0 

207-,  70.0, 
41.5,  20.0, 
6.2 

3 

-58.93 

11.01 

0.0292 

230.,  74.0, 
41.5 

loo.,  180., 
66.0,  14.0 

690.,  200., 
69.0,  24.0 

1* 

27.93 

1.987 

0.0931 

-  .. 

126.,  120., 
91.0,  73*0, 
40.0,  24.0 

-- 

5 

0.0 

2.439 

0.0782 

45-5,  33-0 

70.0,  46.0, 
23.5,  22.0, 

8.9 

-- 
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Qw  cos  8  and  1(6) 


6,  DEGREES 


Figure  10.  Experimental  [l( ©) ]  and  Normalized  Theoretical 
[QnCosG]  Angular  Distribution  Data  for  T  =  30%  L/rm  =  2.0, 

with  Nitrogen. 
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NORMALIZED  Q*  Cos©  AND  1(6) 


Figure  11.  Normalized  Experimental  [ 1 ( 0 ) ]  and  Theoretical 
[QnCos0]  Angular  Distribution  Data  ior  Nitrous  Oxide  through 
Orifice  2  (T  =  8.500%  L/rrn  =  10. 08). 
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NORMALIZED  Qn  Cos©  AND  1(9) 


Figure  12.  Normalized  Experimental  Cl ( 6) ]  and  Theoretical 
[qnCos9]  Angular  Distribution  Data  for  Nitrogen  through 
Orifice  J  (T  -  -58.93°,  L/rm  =  ll.Ol). 


NORMALIZED  Qn  Cos  9  AND  I  (6) 
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Figure  13*  Normalize^  Experimental  [l( 0) 1  and  Theoretical 
[QnCos6]  Angular  Distribution  Data  for  Nitrous  Oxide  through 
Orifice  3  (T  =  -58.93%  L/rm  =  11. Cl). 


NORMALIZED  On  Cos©  AND  I  (6) 


©,  deg 


Figure  lL*  Normalized  Experimental  [l( 6) ]  anJ  Theoretical 
[QnCos9]  Angular  Distribution  Data  for  Nitrogen  through 
Orifice  1  (l  =  25.  (^ ,  L/rm  =  4.010). 


52 


6=0°,  i.e.,  _Iq  ■  Xq/J.0^  with  representing  the  measured  ion  current  at  G. 

In  Figures  9  and  10  the  experimental  points  are  shown,  and  the  lines 
are  drawn  through  theoretical  values.  In  Figures  11-lb,  the  dashed  lines 
represent  smoothed  curves  drawn  through  the  experimental  points  for  both 
positive  and  negative  angles;  the  scatter  of  points  about  the  smoothed  curve 
In  Figure  8  is  t yp i c a 1 .  i.e»,  not  the  minimal.  The  solid  line  in  each 
figure  is  drawn  through  theoretical  values  which  are  discussed  in  the  follow¬ 
ing  section. 

C.  DISCUSSION 


1.  Remarks  on  Theory :  Implicitly  in  Part  I1^  and  CES 13  and  explicitly 
in  Reference  lb  the  theoretical  angular  distribution  is  expressed  as 

dNg( L)  ■«  2rruo  Qn  sin6  CosG  d© ;  (9) 

dNQ(L)  is  the  number  of  molecules  which  effuse  per  second  into  the 
solid  angle  2ti  Sj.nQ  dG,  po  Is  the  molecular  flux  incident  on  the 
entrance  of  the  orifice,  and  _Qn  (n  «*  1,2,5)  is  the  complicated  function 
of  T,  L/xo,  and  G  which  describes  deviations  from  the  "ideal"  cosine 
law  distribution.  We  now  define 


Qv  °  Qn/n  (v  »  n  =  1,2,5),  ( 10) 

insert  Xo  explicitly  (xo  and  i  are  normalized  to  _r0  =  1  In  the  theo¬ 
retical  analysis :l2- 34 ) ,  and  write  in  terms  o£  a  generalized  incre¬ 
mental  solid  angle  dju: 


dNe(L)  =■  ponroc  Qv  Cos©  ( 2n  SinG  d@) 
d^oj(  L )  =  ponr02  Qv  Cos0  du;. 


(11) 


The  transmission  probability  W  may  be  expressed  as  (number  of  molsculet 


r»vf  r* ^  f(  mmiKor  onf  no  ©ri  f  ^  r* ^ ^ 

©r  with  equation  (li) 

W  =  ( l/ponrr^)  J  ponr02  Qv  Cos!  <lu 

( a) 

(12 

■  J  Qv  CosG  dm 

(b) 

“  TqZ  Cos0  (2ti  Sine  de) . 

( c) 

F.quation  (12c)  may  be  compared  with  equation  (l),  if  equation  (10)  is 
noted. 
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We  assume  Chet  Che  aperture  of  the  molecular  beam  detector 
subtends  st  the  effusion  orifice  e  solid  ancle  (q^)  sufficiently  small 
that  £v  Co sS  may  be  considered  constant  over  O-j-  For  two  angles  0 
and  we  may  now  write 


N^(l.)  y ®  nanro*  Qv '  CosO  '  Od 

and 


(a) 

U3) 


H'/  h)  ^  *  ponto"  Qv  Cos0  flj. 

The  rstlc  of  equation  ( 1.3b)  to  (13a)  is,  if  G '  is  taken  ««  0s, 


'dQ(L)c 

Ni/hTo  th 


Qv  Cos 6 

ToTT 


(i) 


(14) 


Equation  (  lH)  gives  the  theoretical  value  of  the  ratio  of  (number 
entering  detector  at  0)  to  (number  entering  M  6  “  C5). 

The  basic  assumption  concerning  the  operation  of  the  beam 
detector  is  that  the  measured  ion  current  _I_+  is  proport ionel  to  the 
number  of  neutral  molecules  entering  the  detector,  or 


"  1%  a  v  ( J5) 

exp 

Wc  r.ow  define  the.  symbol  i  to  mean  nfs  (theoretically)  predicted  to 
be  equal  to",  and  combine  equations  ( 3 4 )  and  (15)  to  obtain 

t0  *  Qv  c°se,  (16) 


4  / 


For  any  orifice,  at  0  =  0°  equation  (26)  becomes 

-  1  *  (Qv)o  •  l  ■»  l;  (IT) 

lo  ( Qv)o 

hence,  both  the  experimental  and  the  theoretical  results  are  self- 
normalizing  to  unity  at  0  ■»  0°.  Equation  (  16)  la  the  basis  for  the 
form  of  the  graphs  in  subsection  B. 

As  indicated  in  equation  (if),  the  plotting  method  suggested  by 
equation  (l6)  and  used  in  Figures  8-1',  forces  Agreement  between 
experimental  and  theoretical  results  at  0  “  0° .  Furthermore,  the 
nature  cf  the  experiment  essentially  forces  agreement  at  0  »  90° • 
Therefore,  whatever  the  actual  nature  of  the  discrepancy  between 


.yi 


experimental  and  theoretical  results,  in  [  lots  of  Jig  atu*-.Qv  CosG/(.,Qv)o 
ye.  G  the  apparent  discrepancy  near  6  =»  0°  is  smell,  and  any  reel 
discrepancies  are  forced  to  appear  in  midrange  of  0  and  are  therefore 
overemphasized.  To  circumvent  this  difficulty  Phipps  and  Adams23  have 
introduced,  and  Wahlbeck  and  Wang27  have  also  used,  a  probability 
density  function  Pg  which  we  now  consider. 

The  experimental  probability  density  function  F*  ie  defined23'37  by 

pe  "  1e+  (2n  st°e  08) 

Obviously,  if  normalized  ion  currents  (ci«  equation  lb)  are  substituted 
for  ,  Pg  ia  unchanged: 

p0  “  1  o/J0TT/s  IG  (2nSinG  d°)  ( 19) 

We  alnc  note  that  since  .I^q  “  1.0, 

i-o  -  v£/a  I6  ^2tt  Sine  de)  (20) 

and 


iepo- 


(21) 


It  is  apparent  from  equation  (2l)  that  «  plot  of  P*  vjs.  p  will 

X 

differ  from  a  plot  of  I.q  vs.  6  onl v  by  the  factor  P(>  For  a  given  set 
of  Qo,0]  data  Po  is  obviously  iixed,  but  there  i  no  requirement  on 
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the  constancy  of  P0  from  run  to  run;  hence,  is  not  self-normali zed 
and  can  reflect  discrepancies  between  experimental  and  theoretical 
data  at  Q  =  0°  as  well  as  at  other  values  of  0- 

The  interpretation  of  Pg  as  a  probability  density  function  follows 
immediately  from  its  definition  (equation  (l8)  or  (19)):  Pq  dun  is 
the  probability  that  an  effusing  molecule  will  transpierce  duj  at  G, 
or  stated  differently,  Pa  is  the  fraction  (cf  cf  fusing  molecules) 


which  ilow,  per  steradian,  at  Q.  With  this  interpretation  of  P 


0 


vif.h  equation  ( 11 )  an  expression  fur  ihe  corresponding  theoretical 


quantity  P.  may  be  obtained: 

pC  =  dHm  (L)  « 
9  dui-Kd)' 


LLonl'o^  Qt  C.OS0  tlj. 

- ” "  -  »«»»■—  ■—  ‘ig— ~ 

UmVoTfio^  W 


or> 


pj  -  (Qv  CO6  0)/W 


(22) 
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The  quantity  Pg  also  is  not  sel i -normalized,  but  may  be  converted 
to  a  self-normalized  form  by  multiplying  by  W/(j)v,)o;  the  result  is 
then  to  be  compared  with  equation  ( 1 6) : 

W  pJ/(Qv)o  -  (QvCosq)/(Qv)0  i  IQ.  (2J) 

From  equetiona  (23)  and  (21)  one  obtains 

I'o  *  (  Qv)o/w  •  (24) 

From  the  above  observations  we  draw  the  following  conclusions: 

(1)  A  plot  of  Pg  vs.  6  reveals  no  information  not  already 

x 

provided  by  a  plot  of  Jig  vs.  0  end  by  the  value  of  P0 
(.cf.  equation  2l) 

(2)  Tf  one  attempts  to  normalize  P*  via  equation  (23)  and  compare 

the  results  with  lQ  (e.g.,  6i  «*  1q  ~  PqW/(Qv)o)j  obtains 

nothing  new,  since 

6 1  -  Ifi  -  l’^/(Qv)o  -  Ie  {1  -  P5w/(Qv)o)  (25) 

and  the  bracketed  portion  of  the  equation  is  a  constant 
for  a  given  run.  If,  on  the  other  hand,  one  attempts 
to  compare  _Ig  directly  with  WFq/(}iv)o  (equation  23)  by 
defining 

6a  =■  [Ie  -  WPj/(Qv)ol/(WPe/(Qv)o], 


(3) 


one  finds  that 

6a  -  [l0(Qv)o/wS]  -  1,  (2b) 

and,  by  virtue  of  equations  (IT)  and  (22),  ta  fa  inherently 
normalized  to  zero  at  0  *»  0° .  Furthermore,  a  plot  of  S2  vs 
0  is,  except  for  a  shift  in  the  zero  point  of  the  ordinate 
caused  by  non-identity  oi  Pq  and  V  2- ,) o/H  (cf.  equation  24), 
idei'iica"  to  the.  more  useful  form  now  to  be  described. 

The  most  informative  scheme  for  comparing  experimental  and 
theoretical  results  is  that  obtained  by  defining,  and  plot¬ 
ting  vs.  e,  the  quantities  A  and  J(a): 


0  ‘  e  . 

-  (rX)  '  1- 


(27) 
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A  plot  of  J[(a)  vs.  9  provides  a  direct,  non-normal i zed 
indication  of  the  discrepancy  between  experimental  and 
theoretical  results,  expressed  as  a  fraction  of  the  theo¬ 
retics  value.  The  only  restriction  on  the  value  of  £(^) 
is  that  Pg0  ™  0  and  Ps0  (i.e.,  J.eo)  is  adjusted  to  zero 
by  subtracting  background  ion  current;  therefore,  i(A)eo*1  0. 
The  one  difficulty  presented  by  the  use  of  _f(  a)  is  that 
at  9  >  -70°  J?q  and.  P^  may  be  rather  small  (e,g„  0.02  of 
the  value  at  9  =  0°);  the  difference  ^  may  be  small  com¬ 
pared  to  experimental  error  in  obtaining  Jg,  but  £(a)  can 
be  quite  large  (e.g.,  0.20  -  0.40)  because  Pq  is  small. 

In  other  words,  the  fractional  errors  are  magnified  as  Py 
approaches  zero.  For  this  reason,  one  might  prefer  to 
plot  £  vs,  9;  both  types  of  plots  are  illustrated  subse¬ 
quently. 

2.  Experimental  Results :  Figures  1S-20  present  various  portions  of 
out  results  in  the  form:  (P^/Pq)  -  1  vs.  0,  i.e.  f_( a )  vs.  9,  except 
for  Figure  16  which  is  a  replot  of  Figure  15  in  the  form  (P*  -  P^)  vs. 
9,  i.e.,  A  vs.  0.  Comparison  of  Figures  IS  and  10  illustrate  the 
point  made  in  the  preceding  section:  the  large  excursions  of  J:(a) 
at  0  >  70°  do  not  realistically  reflect  the  discrepancy  between  P* 
and  Pq.  At  0  <  6o°  plots  of  a  US-  Q  and  of  i_{^)  vs.  9  are  very  similar 
except  for  scaling. 

Figure  21  is  a  replot  in  the  form  [( Ig( Q^)0/WPg)  -  1]  _vs.  9, 
i.e.  vs.  0,  (cf.  equation  26)  of  Run  125  in  Figure  15,  and  is 
included  to  illustrate  two  points:  ( l)  a  plot  of  j,2  vs.  9  is  nor¬ 
malized  to  zero  at  q  =  0°,  and  (2)  except  for  point  ( l)  and  scaling, 
the  shape  of  the  curve  is  the  same  as  that  of  _f(^)  vs.  0  in  Figure  15. 

Before  v* s  sxsnil tiv  tilic  v,  r { ou s  curves  det"!  1.  it  fcc  useful 
to  note  again  the  geometry  cf  the  conical  orifice  and  the  two  important 
angles  T  and  t)*;  these  are  illustrated  in  Figure  22  for  the  diverg¬ 
ing  orifice  (I  >  0)  .  For  0  r£  0  ^  1,  the  detector  may  receive  moleculeB 
di  i  ectiv  flow  of.  the  circular  entrance  to  tin:  o>  lflce  and  from  all 

elements  oi  the  orifice  wall;  this  coj responds  “  to  Range  1  and  n  «  1 
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Figure  ;;2.  The  Critical  Angles  and  the  Angular  Ranges  for 
a  Conical  Orifice 


(c_f.  equa'.icm  9  and  10) ,  For  T  <  9  <  6*,  the  detectoi  may  receive 
molecules  directly  from  elements  oi  the  orifice  wall  over  the  entire 
length  of  the  orifice  and  from  part  of  the  entrance  to  the  orifice, 
but  a  portion  of  the  wall,  and  of  the  entrance,  are  shielded  irom 
the  detector  by  the  outer  rim  of  the  orifice;  this  case  corresponds 
to  Range  II13  and  n  ■  2.  Finally,  for  0*  £  6  £  tt/2  the  detector 
receives  molecules  only  from  part  of  the  orifice  wall;  molecules  which 
traverse  the  orifice  entrance  cannot  proceed  dl rec t ly  to  the  detector; 
this  is  Range  III13,  and  n  ■  3-  For  the  converging  orifice  (T  <  o), 
the  above  statements  are  also  valid  ( l)  if  T  is  replaced  by  |  T  J ,  and 
(2)  If  the  description  of  Range  I  is  modified  to  read;  "receive 
molecules  directly  only  from  „  portion  of  the  orifice  entrance". 

The  significance  of  T  and  0*  with  regard  to  angular  distribution 
may  be  summarized:  For  9  <  | T | ,  j^ll  elements  of  the  orifice  wall 
(T  >  0),  or  jio  elements  (T  <  0),  contribute  to  the  flux  at  the  detector 
aperture,  As  9  increases  in  the  range  J T j  s  9  <  0*,  a  decreasing 
portion  of  both  the  orifice  entrance  and  the  orifice  wall  contribute 
to  the  detected  flux.  For  9  >  9*,  the  orifice  entrance  makes  no 
contribution,  all  molecules  reach  the  detector  from  a  portion  of  the 
orifice  walls,  end  this  portion  decreases  to  0  as  0  -«  tt/2. 

We  shall  now  examine  the  features  of  the  curves  in  Figures  15*20, 
beginning  with  Figures  20  and  15.  The  agreement  between  experimental 
and  theoretical  values  for  Pfi  for  the  k20  series  (Figure  20)  is  very 
good  except  in  the  range  ^0“  i  0  i  60°  where  there  is  a  maximum  in 
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each  curve.  Similarly,  in  the  120  series  the  discrepancy  A  **  Z0  -  Pq 
for  the  lower  pressure  runs  ( 121  and  12j>)  is  zero  within  2  to 
l.e.,  very  nearly  within  experimental  error,  for  0  <  T;  however, 
there  is  again  a  maximum  in  the  curves  at  e  ■  1+0°.  At  higher  pressures 
(Runs  125  and  127)  _f(  a)  is  "0.05  to  -0.10  at  Q  <  T,  the  maximum  in 
each  curve  is  higher  and  is  shifted  to  larger  9.  We  now  note  the 
correlation,  in  Figures  15,  1 6,  and  20,  of  the  maxima  with  the  angles 
T  and  9*,  and  that  the  magnitude  of  both  the  maxima  and  the  discrep¬ 
ancy  1(a)  increases  with  increasing  pressure  and  with  increasing 
(L/r). 

kC 


Support  1  ('!  this  last  gene  i  al  i  zat  ion  was  sought  in  Wang's  data 
for  flic  effusion  of  cesium  chloride  !  hi  ougli  cylindrical  (T  »  0.0°) 
copper  and  nickel  orifices;  data  required  tor  the  ^  _vs.  a  plots  In 

L’|rt.,rf.o  OC  ....U  Ok  1,..,-,,.  „  1  ....  1  i  r.  A  l...  ....  t  , .  n  / 1  1  J  «  « 

»  *-  CjVj  A.  L, ...  1.^1  Qiiu  •_  t  l»uvu  llvtii  lo  i.'}  Uo  1  k  fill  1*0  }  *-  (J  »yj  uol-a 

tabulated  by  Wang1'’ .  Maxima  in  the  various  curves  at  10-15°  and  minima 
at  00-70°  are  immediately  obvious  (but  note  that  the  ^  scale  is  expanded 
by  a  facLor  of  five  compared  with  Figures  15”2l).  The  discrepancy  ^ 
does  indeed  appear  to  decrease  with  increasing  pressure,  but  this  assertion 
must  be  qualified:  P.un  75  In  Figure  2't  is  "out  of  order"  for  no  apparent 
reason;  in  Figure  25  Runs  21  and  2u  were  made  with  one  experimental 
configuration.  Runs  2f,  2fi,  and  55,  with  another;  within  each  group 
A  decreases  with  increasing  pressure;  the  reason  for  the  discrepancy 
between  the  groups  is  not  apparent.  The  shift  of  the  maximum  to  higher 
with  increasing  pressure  is  quite  evident  in  Run  91,  Figure  2^4  ,  but  is 
not  clearly  exhibited  in  Figure  25,  perhaps  because  data  were  not  taken 
at  sufficiently  high  pressures 

It  might  appear  that  Wang's  data,  then, provide  support  for,  if  not 
conf ini'.ation  of,  our  generalization  about  the  variation  of  4  with  (h/_r) 
and  with  pressure.  However,  wc  were  surprised  that  for  low  pressures 
the  maxima  in  Figures  25  and  2lt  occur  in  the  same  angular  range  (10-15°) 
despite  the  difference  in  (L /_r)  for  the  two  orifices.  Upon  investigation, 
we  find  that  in  Wang's  apparatus  v  the  angle  from  the  orifice  to  the  annulus 
which  surrounds  the  baffle  plate  in  the  front  oven  is  in  the  range  10  to 
15°-  It  is  therefore  not  clear  whether  the  maxima  in  the  curves  from 
Wang's  data  arise  from  gas  llow-orificc  phenomena  or  from,  t lie  baffle- 
pnnulus  acting  as  a  (relatively)  concentrated  source  of  molecules. 

An  obvious  question,  then,  is  whether  the  maxima  in  Figures  15,  Itj 
and  20  arise  from  c  cell  effect  rather  than  an  orifice  effect.  The 
cylindrical  interior  of  our  simulated  Knudsen  Cell  Is  1.00"  in  diameter 
and  l.lo"  high;  the  bottom  of  the  cell  is  completely  open  to  the  flow 
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t]i«  ir.t  erior  wsll  of  the  cell  Is  rnscro- 


scopically  smooth.  Hence,  the  only  discontinuity  in  the  emitting  surface 
"seen"  by  the.  detector  through  the  effusion  orifice  is  at  the  "bottom 
corner"  of  the  cell;  the  angle  from  the  orifice  to  the  "bottom  corner" 


7 


is  -25° •  For  o  greater  than  —29°  molecules  can  travel  directly  i r <xu 
the  cell  wall  through  the  orifice  (if  0*  >  23°)  to  t lie  detector;  for  0 
less  than  ~£5° >  the  detector  cannot  ’’flee"  the  cell  wall  and  receives  no 
molecules  directly  from  it.  Wc  now  examine  Figures  19  and  20  for  evidence 
ci  a  cell  effect  an  find  that  in  Figure  10  jf  is  so  near  25°  that  one 
cannot  hope  to  distinguish  cell  effects  from  orifice  effects;  in  Figure 
20  the  situation  is  similar.  However,  the  data  for  Figures  13  and  19 
were  obtained  for  a  converging  orifice  with  T  -  ~5S.9°>  in  this  case  the 
orifice  presents  no  obstruction  to  molecules  traveling  along  a  6  *  25° 
trajectory  from  the  region  of  the  surface  discontinuity  to  the  detector. 
Figures  17  and  18  provide  no  evidence  (i.e.,  no  maxima,  no  minima,  no 
breaks  in  the  curves)  for  a  cell  effect  at  or  near  Q  -  25°,  but  there 
is  a  definite  change  (more  pronounced  at  higher  pressures)  in  the  curvature 
of  the  various  curves  in  the  narrow  range  |t|  s  g  s  0*  .  We  therefore 
conclude  that  the  deviations  from  theoretical  behavior,  exhibited  by  the 
curves  in  the  various  figures,  are  ori flee  effects. 

We  return  our  attention  to  Figures  10  and  20  and  consider  the  two 
mechanisms  frequently  cited  as  causes  of  deviation  from  theoretical 
behavior:  specular  reflection  from  the  orifice  wall  and  surface  diffusion 
along  the  orifice  wall0.  One  would  expect  any  enhancement  of  intensity 
by  specular  reflection  to  occur  at  angles  0  T,  Enhancement  at  angles 
0  >  T  requires  apeculer  reflection  through  angles  greater  than  2T,  and 
the  maxima  of  Figures  16  and  20  would  require  preferential  enhancement 
by  specular  reflection  through  angles  >  50°  with  no  concurrent  depletion 
at  smaller  anglea  (_c£.  especially  Figure  20) .  Available  experimental 
date39  on  angular  distribution  of  molecules  reflected  from  surfeces  do 
not  support  these  rather  stringent  requirements. 

If  aurface  diffusion  occurs  along  the  orifice  wall,  the  concentration 
of  molecules  resident  within  any  incremental  area  on  the  wall  ia  expected 
to  be  greater  than  if  aurface  diffusion  did  not  occur;  consequently,  the 
flux  from  the  orifice  wall  to  the  detector  is  expected  to  be  greater. 
Therefore,  one  might  expect  any  discrepancy  between  theoretical  and 
experimental  data  for  angular  distributions,  which  results  from  surface 
diffusion,  to  be  moat  prominent  in  the  range(s)  of  0  in  which  the  flux 
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from  the  orifice  wall  area  makes  l he  maximum  tclexfvc  contribution  to  the 
total  flux  into  the  detector. 

The  area  of  the  projection  of  the  orifice  wa 1 1  onto  a  plane  per* 
pendicular  to  the  orifice-detector  line  defined  hy  p  may  be  calculated 
from  the  following  equations  (see  Figure  1  for  symbols): 


Ag  "  ni'k"  Cos  0 

Aq  ■  nro,!  Cos  c  (fO) 

Ap  »  A(  IIq  1  Dp,)  Cos  0 

("A( D0 1  Dg)  is  the  area  of  overlap  of  the  circle  nto*  with  the  pro 
jection  at  angle  p  onto  the  plane  of  that  circle  of  the  circle  ^r^2; 
derivation  of  this  quantity  has  been  described  in  detail 12' 14.  ] 

if  e  v»  e*,  aw  -  al 

If  |T|<  0  <  0*.  Aw  =  Al  -  Ae  (Oi) 

If  6  s  |l'|  and  T  >  0.0,  Aw  =  Ag  -  A0 
If  o  j T |  and  T  <  0.0,  Aw  “0.0 


Th.e  results  of  these  calculations  for  the  criliecs  for  which  u..gv.l„x 
distribution  data  were  obtained  arc  shown  in  Figures  25-2 0  wherein, 
for  convenience  in  plotting,  the  areas  have  been  normalized  to  Ag(  o“0. 0)**1. 0. 
In  the  ranges  p  ;>  p«  and  p  -±  T,  the  variation  of  Ay  with  9  is  given  by 
Cos  0.  However,  in  the  range  T  •;  q  s  c*,  Aw  exhibits  a  rather  different 
behavior  which  produces  a  minimum  (ox  virtual  minimum)  in  Aw  at  9  «  T 
and  a  maximum  (or  virtual  maximum)  at  9  “  9 . 


As  0  increases  over  the  range  T  •:  9  ^  9*  lor  a  diverging  (T  >  0.0) 
orifice,  the  entrance  of  the  orifice  is  eclipsed  by  the  enter  rim  of  the 
exit,  and  therefore  the  fracti onal  contribution  from  the  wall  to  the 
total  flux  to  the  detector  increases  to  unity  at  9  =  p-H.  In  the  range 
(9  S  0W)  over  which  the  fractional  contribution  iron-,  tlie  wall  is  unity, 
the  projected  wall  area  "seen"  by  t lie  djtector  Is  a  maximum  at  q  =  (d*-. 
Therefore,  _i_f  there  were  an  increase  in  the  flux  from  the  orifice  wall 
above  that  predicted  by  our  extension  of  Clnusing's  Model,  one  should 
expect  to  see  maxima  in  ^  ve.  9  plots  <  t  or  near  0  =  g-* . 

This  is  precisely  what  is  observed  in  Figures  lp,  lo,  1Y,  and  20,  and 
in  lieu  of  any  acceptable  alternative,  wc  tentatively  conclude  that  the 
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increased  fluK  from  the  orifice  wall  results  from  surface  diffusion. 

We  have  not  yet  made  a  more  detailed  analysis  of  a  model  including  surface 
diffusion;  Ruth  and  Hirth29  have  reported  the  corresponding  analysis  for 
cylindrical  orifices. 

Other  features  oL  the  j[(a)  vs.  q  plots,  e,g.  variations  with  pressure 
ar.d  nature  of  the  gas,  are  being  analyzed  and  will  be  reported  in  the  open 
literature. 
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SECTION  IV 

THE  MULT I CELL  TECHNIQUE  FOR  EXPERIMENTAL 
DETERMINATION  OF  TRANSMISSION  PROBABILITIES  FOR 
MOLECULAR  FLOW  THROUGH  CONICAL  ORIFICES 

The  Knudscn  effusion  technique1"3  Las  been  widely  used  to  obtain 
vapor  pressure  data  in  the  range  10'3-lL)'0  atm.  The  technique  has  several 
variants3,  but  basically  the  experimentally  measured  variable  is  the 
isothermal  rate  o£  flew  of  vapor  through  an  orifice,  e.g.,  1  mm  in  diameter 
and  0.5  nvn  long.  The  equation  which  relates  the  pressure  of  the  gas  to 
the  flux  through  the  orifice  contains  a  factor  which  is  often  called  the 
Clausing  factor,  or  more  properly,  the  transmission  probability,  of  the 
orifice,  l.e,,  the  probability  that  a  molecule  which  has  entered  one  end 
of  the  orifice  will  exit  from  the  opposite  end.  Numerical  values  for  the 
transmission  probability  are  available  from  theoretical  analyses5  B> 13  14 
of  molecular  flow  of  gas  through  orifices  with  various  geometries.  Although 
these  theoretical  transmission  probabilities  are  widely2'3  used  to  collect 
the  observed  rate  of  flow  for  the  effect  of  finite  orifice  length,  reliable 
experimental  data  with  which  to  compare  the  theoretical  transmission  prob¬ 
abilities  are  scarce3*31/33. 

The  Multicell  technique,  which  we  describe  below,  provides  a  simple 
method  for  obtaining  experimental  values  for  transmission  probabilities 
and  also  offers  a  means  for  critically  testing  various  assumptions  which 
have  been  made9'34”36  concerning  molecular  flow  within  a  Knudscn  cell. 
Application  of  the  technique  to  the  determination  of  experimental  transmission 
probabilities  is  discussed  in  this  section. 

A.  THE  MULTICELL  APPROACH 

In  a  Multicell  effusion  experiment  a  number  (e.g.,  <+  to  10 )  of  ef¬ 
fusion  cells2’3  which  are  as  nearly  identical  as  possible  excep  t  in  the 
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variable  under  study  are  subjected  simultaneously  to  a  given  environment 
for  a  given  time.  Comparison  of  results  (c.g.,  weight  lost  from  each  cell) 
for  the  Several  cells  _in  a  given  experiment  permits  isolation  of  the 
effect  of  the  variable  of  interest  since  all  other  variables  are  constant, 
or  vary  simultaneously  and  identically,  for  all  cells. 

In  the  application  of  this  approach  to  the  determination  of  experimental 
transmission  probabilities,  Knudsen  c«_lls  are  fabricated  to  be  essentially 
identical  except  that  each  orifice  is  different  from  the  others;  each  cell 
is  loaded  with  an  equal  amount  of  a  given  sample  (e.g. ,  cadmium  or  mercury), 
and  then  heated  simultaneously  in  vacuum  for  a  given  time  T  and  at  an 
appropriate  temperature  Q  which  sre  identical  for  all  cells.  The  weight 
lost  from  each  cell  during  the  experiment  is  determined  by  weighing  before 
and  after  the  heating  period.  The  weight  d^  lost  from  the  _ith  cell  In 
time  ,dt  Is  given  bynT 

dgl  o  aiK1[P(M/2nRo)2y/;3dt]  (JO) 

in  which  JP  is  the  pressure  in  the  cell  at  temperature  0,  M  is  the  molecular 
weight  of  the  effusing  molecules,  is  the  experimental  transmission 
probability  and  the  cross-sectional  area  for  the  jLth  orifice.  For 
the  experimental  conditions  described  above,  P,  0,  and  M  are,  at  each 
instant,  the  same  for  all  m  cells;  therefore,  the  integral  of  the  bracketed 
portion  of  equation  (30)  over  the  interval  0  to  T  is  the  same  for  each  cell. 
Hence,  if  is  the  total  weight  lost  In  time  T,  we  write 

.ft=0  dgl  =  gi  "  eiKfi 

with  1  representing  the  integral  which  is  constant  for  the  several  cells 
in  a  given  run,  but  which  varies  from  run  to  run.  Since  the  various 
and  are  directly  measurable  quantities,  it  is  convenient  to  define 
■  JJj/ili  and  then  to  write 

Ji  -  KtI.  (51) 

Clausing's  theory  predicts5"8?  13  that,  for  conical  (-n/2  <  T  <  tt/2) 
and  cylindrical  (T  -  0)  orifices,  the  variation  in  a  given  run  of  the  values 
j  for  the  various  cells  depends  only  on  the  ratio  of  length  L  to  diameter 
D  for  the  orifices  and  on  T  the  semi-apex  angle  of  the  cone  of  which  the 
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orifice  is  8  truncated  section.  .  Iietsforo,  if  tl'c  values  Jp  are.  plotted 
against  some  suitable  function  of  (b/p)  and  T,  designated  by  J.’(.b  'P.T) ,  end 
extrapolated  to  Tv  0,  T),  one  ruay  write  from  Equation  (51),  Jo  “  BcJ,  which 
the  subscript  designates  t lrc  condition;  ]•■'  1)  =  O'.  But  (b'D  =  0)  defines5  M>  13 
tic  itidcale  orifice  for  which  the  transmission  probability  b0  is  unity 
legardless  of  the  value  of  T  (which  in  fact  is  mesningless  when  L/D  =  0); 
therefore, 

Jo  -  I  (32) 

Hence,  the  intercept  at  K(0,  T)  is  just  the  weight  loss  per  v - i t  area  J0 
to  be  expected,  under  the  conditions  of  the  run,  from  an  Ideal  orifice. 

With  this  experimental  value  for  Jl ,  the  experimental  transmission  probability 
for  each  of  the  cells  is  calculated  irom  Equations  (jl)  and  (}2): 

Kp  -  Jt/J0-  (33) 

The  most  suitable  choice  for  F(L/B,T;  appears  to  he  W,  the  theoretical 
transmission  probability  s"b> 13  calculated  from  the  orifice  dimensions.  This 
choice  results  in  a  conveniently  bounded  plot  with  abscissa  values  ranging 
from  zero  to  one;  more  important,  if  the  theoretical  transmission  probabilities 
predicted  by  Clausing's  theory5  and  its  extensions6-0?  13  are  correct:,  a  plot 
of  versus  F(L/D,t)j  «  Wj_  should  give  a  straight  line  which  may  be  precisely 
extrapolated  to  =  1  (i.e.,  to  L/D  -  0)  to  obtain  the  irtarcept  J.0-  The 
slope  of  the  straight  line  resulting  from  such  a  plot  should  also  equal  Jo- 

B.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

The  results  reported  here  were  obtained  with  a  typical  Multicell 
configuration:  various  sets  of  eight  nichel-plated,  steel  Knudsen  cells, 

machined  to  be  "identical"  except  for  the  lengths  of  the  cylindrical  orifices 
and/or  the  semi-apex  angles  T  of  the  conical  orifices,  were  placed  into 
closely-fitting,  symmetrically-arranged  cavities  in  an  aluminum  (or  copper) 
block  which,  during  a  run,  sat  inside  a  stainless  steel  vacuum  chamber. 

During  the  heating  period,  the  entire  vacuum  chamber  and  several  inches  of 
connecting  stainless  steel  vacuum  line  were  inside  a  circulating-hot-air, 
constant- temperature  oven.  The  cells  were  charged  with  equal  amounts  of 
high-purlty  cadmium  metal. 


The  dimensions  of  the  orifices  and  the  detailed  results  of  the  various 
runs  are  available  elsewhere30* 40 j  the  results  of  a  typical  run  are  presented 
in  Table  3  and  in  Figure  29  in  which  the  least-squares  straight  line  (j^C  “ 
mW^  +  J>)  for  the  data  points  is  also  given. 

For  each  run,  the  experimental  transmission  probabilities  were  calculated 
in  two  ways:  (l)  the  ordinate  value  of  eacli  experimental  point  (i.e.,  each 
was  divided  by  _J0  (cf.  equation  33)  tp  obtain  the  values  K^;  (2)  the 
cori esponding  ordinate  value  (j^c)  jnn  t he  least-squares  line  was  divided  by 
J0  to  obtain  the  values  KjC.  In  Table  3,  the  transmission  probabilities  cal¬ 
culated  in  the  two  ways  from  the  data  of  Run  7  are  compared  with  the  theoretical 
transmission  probabilities. 

For  each  orifice  the  transmission  probability  calculated  from  the  results 
of  several  runs  will  be  more  reliable  than  a  value  obtained  from  any  single, 
arbitrarily-chosen  run.  "Average"  transmission  probabilities  for  the 
various  sets  of  orifices,  and  the  transmission  probabilities  _K^  calculated 
from  the  individual  runs  are  given  in  Tables  4-6.  Two  values  are  given  for 
each  orifice:  the  upper  of  each  pair  of  numbers  is  the  transmission  probability 
calculated  from  the  experimental  poin*  ;■  (cf.  K^p,  Table  3),  and  the  lower 
number  is  the  value  calculated  from  the  least-squares  line  for  the  points  of 
a  given  run  (.cf.  K^C,  Table  3 )  • 

As  shown  in  the  last  column  of  Tables  4-6,  the  discrepancies  between  the 
experimental  values  for  transmission  probabilities  (  foi  cadmium  vapor  passing 
through  nickel-plated,  cylindrical  orifices)  and  the  theoretical  values  lie 
in  the  range  ±3^  (with  two  notable  exceptions  iri  Table  6,  for  which  we  have 
no  explanation),  if  the  experimental  points  are  considered.  If  the  results 
from  the  least-squares  lines  are  used,  the  discrepancies  are  appreciably 
smaller.  Similar  agreement  for  a  very  limited  number  of  "short"  orifices 
lias  been  obtained  by  McKinley  and  Vance33  and  by  Carlson,  Gilles  and  Thorn37. 

In  general,  the  experimental  values  are  less  than  the  theoretical 
Wj,  and  there  is  a  alight  trend  for  the  discrepancy  to  become  larger  as  L/D 
increases.  These  observations  are  exactly  opposite  to  what  would  be  expected 
if  surface  diffusion0* 36  through  the  orifice  contributed  significantly  to  the 
total  efflux. 

The  agreement  between  theoretical  predictions  and  the  experimental 
results  reported  here  (especially  Table  4)  is  appreciably  better  than  the 
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TABLK  3 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 
'.'■/•..OMISSION  PROBABILITIES:  RUN  7 


Orifice,, 


L/D 

Mi 

kP 

0.195 

O.837 

42.36 

0.843 

0.290 

0.777 

38.16 

0-759 

0.414 

0.711 

35.62 

0.70  9 

0.458 

0.695 

52.94 

0.656 

0.521 

0.664 

33.83 

0.673 

0.522 

0.663 

31.91 

0.635 

0.667 

0.609 

28.85 

0. 574 

0.800 

0.567 

28.55 

O.568 

1  oof  kD-w) 
V’ 

Kc 

ioo(kc-w) 

w 

+0.7 

O.83I 

-0.8 

-8.3 

O.767 

-1.3 

-0.3 

0.699 

-1.8 

-6.0 

0.682 

-1.8 

+1.4 

0.649 

-2.2 

-5.4 

0.649 

-2.2 

-6.1 

0.592 

-2.9 

+0.2 

0. 549 

-3.4 

TABLE  7 


EXI'ER I MENTAL  TRANSMISSION  PROBABILITIES 
FOR  FICFT  CTLIMIM'.’C.'.F  ORIFICES  (SET  l) 


Orifice, 

L/D 

w 

K 

K 

ioo(k-w) 

Run  it 

Run  5 

Run  6 

Run  7 

w 

0. 195 

0.8>r 

0.855l|a 

O.8512 

O.8550 

0.87.29 

0.876 

+1.1 

0. 82517  b 

O.8592 

O.8578 

O.8505 

0.837 

-0.7 

0.290 

0.77  r 

— 

0.7207 

0.7576 

0.75513 

0.775 

-7.  ] 

— 

0.775J1 

0.7772 

0.7072 

0.775 

-0.3 

c.717 

0.711 

0.69j2° 

0.7595 

0 . 6888 

0.7088 

0.708 

-0.7 

0.0975b 

0.7173 

0.7116 

0. 69S6 

0.706 

-0.7 

0.  48 

0.695 

0.665t>f 

O.6018 

0.6951 

0.6555 

0.678 

-2.7 

0.68l0b 

O.6986 

0.6500 

0.0027 

O.090 

-0.7 

0.521 

0. 66  it 

0.6310° 

0.7715 

0.6905 

0. 6752 

0.689 

+3.8 

0.6777b 

0. 6672 

0.6675 

0.6792 

0.657 

-1.1 

0.522 

O.665 

0.6617° 

0.6572 

0.6828 

0.6550 

0.659 

-0.6 

0-6^69 

0.6665 

0. 665 6 

0.6785 

0. 658 

-1. 1 

0.667 

O.609 

0.  S7Q6° 

0.6120 

n  5205 

r\  r.  )  n 

r'  iOo 
v.>. 

-5.3 

0. 5905 b 

O.615O 

O.6096 

0.5921 

0.601 

-1-3 

0.800 

0.567 

0.55^9f 

0.5556 

0. 5701 

0.5681 

0.537 

-1.8 

0. 576>b 

0.5716 

0. 5679 

0.5785 

0.559 

-1.7 

in  =  slope 

58.15 

3-7.57 

80.76 

52.55 

b  =  intercept 

at 

(wt  =  0) 

-1.71 

+O.58 

+9.  20 

-2.  m 

Jrp  intercept 

at 

(wt  ■=  1) 

56.lt! 

57 . 95 

00. 

50. 26 

a.  The  upper  number  of  each  pair  is  (.lie  value  obtained  from  the  experimental 

points . 

b,  The  lower  number  of  each  pair  i3  t  lie  value  obtained  from  the  least 

squares  line. 
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TABLE 


EXPERIMENTAL  TRANSMISSION  PROBABILITIES 
FOR  EIGHT  CYLINDRICAL  ORIFICES  (SET  II) 


Orifice, 

l/r 

w 

Run  2 

Run  3 

0.4462 

0.819 

0.877 a 
0.848b 

0.854 

0.811 

0.66l4 

0.754 

0.735a 

0.761b 

0.752 

0.742 

1.016 

0.669 

0.647a 

0.680b 

0.649 

0.652 

1.291 

0.616 

0.6l4a 

0.603b 

0.594 

0.593 

1.627 

0.563 

0. 590a 
0. 575b 

0-535 

0.541 

1.946 

0.321 

0.536a 

0.535b 

0.484 

0.497 

2.26.5 

0.485 

0.508a 

0.501b 

0.456 

o.46o 

2.573 

0,455 

0.464a 

0.469b 

0.448 

0.428 

m  m  slope 
b  -  intercept 

46.67 

66.15 

fit  (Mi  ■  0) 

J rT  intercept 

+1.24 

-3.10 

at 

(Wt  -  1) 

47.91 

63.05 

1Q0( K-W) 


Run  4 

Run  6 

K 

w 

0.809 

0.797 

O.0l4 

0.802 

O.838 

0.815 

+2.3 

-0.5 

O  O 

ro 

0.751 

0-731 

0.748 

0.740 

-0.8 

-1.9 

0. 593 
0.629 

0.599 

O.638 

0.622 

0.650 

-7.0 

-2.8 

0.541 
0. 570 

0.577 

0.595 

O.582 

0.592 

-5-5 

-3.9 

0.497 

0.510 

0. 520 
0.523 

0. 536 
0.537 

-4.8 

-4.6 

0.482 

0.464 

0.486 

0.477 

0.497 

0.493 

-4.6 

-5-4 

0.419 

0.423 

0.4o4 

0.437 

0.447 

0.455 

-7.6 

-6.2 

0.414 

0.390 

0.446 

0.405 

0.443 

0.423 

-2.6 

-7-0 

84.28 

73-14 

-9-07 

-6.12 

75-21 

67.02 

a.  The  upper  number  of  each  prir  Is  the  value  obtained  from  the  experimental 

points. 

b.  The  lower  number  of  each  pair  is  the  value  obtained  from  the  least 

squares  line. 
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TA1S1.];  u 


EXPERIMENTAL  TRANSMISSION  l’KOHARlLJT  1ES 
FOR  EIGHT  CONICAL  ORIFICES  (SET  III) 


Orifice, 

L/R 

T,  deg 

n 

K 

AVLlil  U 

Run  9 

- 

mn(  k-w) 
w 

1.0S0 

22.  0 

0.8833 

0.837a 

0.650 

0.844 

-4.7 

0.887b 

O.863 

0.886 

»-u.  1 

3.889 

26.8 

0.8848 

0.942a 

0.860 

0.901 

+1.8 

0.880b 

0.884 

O.885 

<0.1 

1.003 

18.7 

0.8665 

0.832a 

0.936 

0.384 

+2.0 

0.867  b 

0.662 

0.864 

-0.3 

0.4218 

0 

0.8252 

0.870a 

0.837 

0.864 

+4.7 

0  827  b 

0.825 

0.826 

+0.1 

4.031 

15.2 

0.7375 

0.620a 

0.581 

O.oOO 

-18.7 

0.742b 

0.736 

0.739 

+0.2 

2.183 

9.63 

O.6898 

0.877a 

0.829 

0.853 

+23.6 

0.695b 

0.688 

0.692 

+0.3 

3.837 

9-15 

0.6187 

C.568a 

0. 595 

0.582 

-6.0 

0.625b 

0.617 

0.621 

+0,3 

2.072 

0 

0.5059 

0.505a 

0.498 

0.502 

-0.8 

0.515b 

0.  504 

0.510 

+0.8 

m  <=>  slope 

37 .91 

35.88 

b  =■  intercept  at 

(Wjoo) 

+0.68 

'0. 17 

Jo3  Intercept  at 

(Wi-l) 

58.65 

35.71 

upper  nuns  be  1  0  i 

- -  1,  „  _  J  - 

vovn  pan 

-  J  ~  ♦  l.  .. 

X  9  t  lie 

V a  lu  e 

—  a  r  — 
vtaamtu 

points, 

b.  The  lower  number  of  each  pair  is  the  value  obtained  from  the  least 
squares  line. 


typical  agreement,  between  measured  vapor  pressures  reported  by  two  laboratories, 
both  of  which  use  supposedly  the  same  Knudseu  effusion  technique,  and  is 
better  than  the  reproducibility  in  vapor  pressure  measurements  often  obtained 
within  a  given  laboratory.  It  would  therefore  appear  that  the  source  of 
these  discrepancies  in  vapor  pressure  measurements  may  lie  in  phenomena  other 
than  flow  of  gas  through  an  orifice. 

Concurrently  wiLh  our  work,  Macur,  Edwards,  and  Wahlbeck41  have  attempted 
to  circuc"'nt  Lhis  problem  ot  poor  reproducibility  in  Knudsen  effusion  measure¬ 
ments  by  using  a  Multicell  technique  in  which  the  emphasis  is  on  the  simul¬ 
taneous  determination  of  several  vapor  pressures  at  the  same  temperature;  the 
avetage  deviation  in  the  pressures  calculated  from  the  several  cells  in  a 
given  run  range  from  1  to  for  indium  and  from  1  to  Qf.  tor  gallium.  The 
agreement  from  run  to  run  (i.e.,  with  change  in  temperature),  as  reflected 
in  the  calculated  values  of  of  vaporization,  is  equally  impressive; 

the  uncertainty  quoted  is  10. 10  Real/ mole. 
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SECTION  V 


THE  MIRER  TECHNIQUE 

The  Miker  technique  for  determination  of  vapor  pressures  and  of  molecular 
weights  of  vapors  at  high  temperatures  \%’as  described  in  detail  in  Part  Ilrl. 

In  this  technique  a  vacuum  microba lance  is  used  to  obtain  both  the  rate 
of  effusion  from  a  suspended  Knudsen  cell  and  the  recoil  force  which  is 
exerted  on  the  cell  ss  a  result  of  effusio,.. 

The  results  reported  in  Part  II"1  demonstrated  the  general  validity 
of  the  technique,  and  also  demonstrated  the  need  for  various  refinements 
if  the  technique  was  to  provide  a  significant  Improvement  over  typical 
Torsion- Knudsen  of fusion-Rc^ui 1  data.  The  various  refinements  will  be 
reported. 

A.  MODIFICATIONS  TO  APPARATUS 

The  early  results21  with  the  Hiker  apparatus  indicated  two  major  de¬ 
ficiencies:  ( l)  When  power  to  the  furnace  was  terminated  in  the  course 
of  measuring  the  recoil  force,  the  furnace  cooled  so  slowly  that  an  un¬ 
desirably  large  correction  for  effusion  during  cooling  was  required. 

(2)  Manual  control  of  the  vacuum  microbalance  was  satisfactory  for  rate 
of  effusion  measurements,  but  not  for  recoil  force  measurements;  automatic 
balance  control  appears  to  be  a  requirement  for  valid  recoil  force  measure¬ 
ments.  Most  of  the  modifications  to  be  described  were  made  in  the  course 
of  overcoming  these  two  deficiencies. 

1.  Furnace  and  Power  Supply :  Rapid  cooling  by  the  furnace  appeared 
to  be  most  easily  achieved  through  the  use  of  elements  with  very  low 
heat-capacity  and  with  high  surface  area.  Tungsten  mesh  furnace 
elements  were  considered  but  rejected  in  favor  of  graph! t c  t ape  for 
the  following  reasons:  (l)  Graphite  tape  retains  its  flexibility 

r." 


alter  in, my  cycles  i  .>  npewu  ing  t  •.•niju-x'.'-- 1  xii't's  ;  (2;  Required  fui  iidtc 
t  emporn  t  uivs  were  not  ,  a  t  i  cast  initial  1  y ,  so  high  t  Ira  t  vapor izat  ion 
Jri  ow  t  he  graphite  el  nient  s  posed  a  p  rob  1 «.  tn ;  (3)  Furnace  fabrication 
is  s  impl  i  1  ied  and  is  readily  aecoaip  1  is  bed  in  our  shop;  (4)  The  resis¬ 
tivity  of  graphite  is  so  high  that  one  may  easily  fabricate  a  furnace 
which  operates  directly  from  an  SCR-eont  rol  led ,  IC'OA,  1 10V,  6011/  line; 
power  at.  this  volt  age- current  level  is  much  more  easily  handled  and 
introduced  into  the  vacuum  system  than  is  the  lower  voltage  -  very 

high  current  levels  required  by  metallic  heat  in;;  elements,  { 

The  furnace  has  a  cylindrical  configuration  and  consists  of  six  13- cm 
lengths  of  1  .27-em-wide  graphite  tope  held  at  top  and  button:  by  machined  | 

graphite  rings  with  graphite  screws.  Electrically,  the  three  lengths 

of  tape  which  form  one  hemicyl  nth  r  are  in  parallel,  as  arc  the  other  | 

throe;  the  two  sets  of  three  are  then  in  series,  connected  electri-  j 

cally  hv  the  bottom  grapni te  ring.  The  resistance  of  the  furnace  is  I 

,V1  .0  ohm  at  2.V'IC,  With  no  radiation  shielding,  a  current  of  ^.50A  rms  5 

produces  a  temperature  of  1300°  0  in  a  ceil  suspended  in  the  ccnte, 

of  the  furnace.  When  power  to  the  furnace  is  terminated,  the  suspended  ■  ' 

cell  cools  initially  at  a  rate  of  ->.25-30  deg  C/sec.  ( 

As  indicated  above,  power  to  the  furnace  controlled  with  a  ,  | 

i 

single  1 10A  si  1 icon- coni rol led- rectifier  the  gate  of  which  is  driven 

4  3 

by  a  "Sil  i  control"  pulse  generating  unit  . 

The  other  units  in  the  control  loop  are  a  t  tings  ten- rhenium  II. er-  5 

mocouplo  which  responds  to  changes  in  furnace  temperature,  a  stable, 
variable  reference  voltage,  an  electronic  null  detector  which  responds 
to  any  error  voltage  between  the  thermocouple  and  the  reference  voltage, 

and  a  Controller  which  accepts  the  output  of  Ihe  null  detector,  provides  • 

proportional,  derivative,  and  integral  control  actions,  and  supplies  an 
appt  opriate  dc  current,  to  drive  the  Si  3.  i  control  unit.  The  thermocouple 
is  used  only  for  control;  temperatures  are  me  ns  tired  with  an  automatic  photo¬ 
electric  "i'lioL  ornj  t  i  <■"  pyrometer. 
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2.  Ant nma tic  Con t ro  1  of  Mi  croho  1  once  :  Several  changes  have  been  made 
in  the  automatic  control  system,  although  the  basic  control  loop 
(Figure  20  in  reference  21)  remains  the  same.  The  magnetic  amplifier 
which  supplied  current  to  the  compensating  coils  has  been  replaced 
with  a  we  1 1- xegula t ed  dc  power  supply  and  an  emitter- fol lower  circuit 
which  is  driven  by  the  output  of  the  Controller.  This  substitution 
eliminated  a  significant  ripple  voltage  from  the  coils. 

The  "Fhotopot"  light  beam  position  sensor**'  has  been  replaced 
with  one  of  later  manufacture,  which  is  more  stable  and  which  can 
he  readily  operated  in  high  vacuum.  The  new  "Photopot"  was  mounted 
ins ide  the  balance  chamber  about  6.3  cm  from  the  mirror  affixed  to  the 
center  of  the  balance  beam.  The  light  beam  originates,  outside  the 
vacuum  system,  in  a  50W  projector  lamp,  is  defined  by  a  0.010-cm  slit, 
reflected  through  a  lens  and  window  into  the  balance  chamber  and  onto 
the  balance  mirror,  and  focused  onto  the  "Photopoi".  The  output  of 
the  projector  lamp  is  monitored  and  controlled  by  a  modification  of 
the  circuit  described  by  Rosencha T 3 .  Power  ior  the  "Phot.opot"  bridge 
circuit  is  now  supplied  by  mercury  batteries. 

The  weight  change  which  can  be  compensated  by  change  in  current 
through  the  balance  coils  (from  0  tc^.lOOmA)  is  28mg.  To  avoid 
opening  the  vacuum  system  to  re-tare  the  balance  after  each  25-30tng 
weight  loss,  a  taring  mechanism  operable  from  outside  the  vacuum 
system  was  installed.  A  bel lows- sea  led  shaft,  driven  outside  the 
vacuum  chamber  by  a  micrometer  screw,  carries  four  -^,20-mg  ring  weights 
formed  from  Alumel  wire  0.020  cm  in  diameter.  As  the  miciomeler  screw 
is  rotated,  the  shaft  is  translated  downward  (upward)  and  the  four 
ring  weights  are  sequentially  placed  on  (removed  from)  a  cross-arm 
attached  to  the  balance  beam. 

3.  M  i  sc  o  1  laiicous  Mod  i  f  i  ca  t  ions:  In  the  Miker  technique  effusing 
vapor  is  directed  downward,  directly  toward  the  window/prism  Lhrough 
which  optical  pyrometer  measurements  are  mode.  A  shutter  protects 
the  window  during  the  interval  between  pyrometer  reading;  even  so, 
the  window  acquires  a  significant  deposit  rather  rapidly.  To  avoid 
o, ening  the  entire  vacuum  system  t  c<  the  atmosphere  while  the  window/ 
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1 

I 

* 

1 

* 

prism  is  removed ,  cleaned,  and  replaced,  very  tbinf  hellows- 
sealed  gate  vulve  was  Installed  in  Hit*  window  mount  and  between  the 

? 

v  indow/j'i  i  s:r.  itself  mrf  the  furnace  cl.imbcr.  The  valve  is  a  modlfi-  , 

3C?  t  4 

cation  of  the  design  reported  by  Sheffield  '  ,  and  has  performed 
sat  i sf act  or t ly . 

The  effect  of  building  and  equipment  vibrations  on  the  micro- 
balance  has  been  minimized  by  mounting  the  vacuum  system  on  a  steel- 
reinforced  concrete  block  (2x2x3  ft.  high)  which  was  cast  in  place 
on  the  basement  floor;  mechanical  vacuum  pumps  in  the  vicinity  have 

been  mounted  on  vibration-absorbing  pads,  etc.  ! 

For  increased  pumping  capacity,  there  has  been  installed  in  the  * 

vacuum  system  a  larger  diffusion  pump  (PMC- 720)  and  a  liquid  nitrogen  j 

trap  (Cryo- sorb) .  j 

4.  Spurious  Recoil  Forces :  In  measurements  of  rate  of  effusion  and  : 

of  recoil  force  made  immediately  after  the  modifications  described  j 

above  were  effected,  the  measured  recoil  force  was  much  larger  (e,g.,  j 

by  a  factor  of  1. 5-2.0)  than  expected  from  the  rate  of  effusion  .  | 

measurements.  A  major  source  of  tire  spurious  portion  of  the  recoil  j 

force  has  been  identified  as  an  interaction  between  the  fields  .  ] 

surrounding  the  furnace  power  leads  and  the  magnet  attached  to  the  ' 

balance  beam.  This  interaction  has  been  minimized,  and  all  but  elim-  j 

lnated,  bv  arranging  the  furnace  power  leads  and  the  furnace  elements  j 

themselves  in  a  configuration  symmetrical  with  respect  to  generation  J 

ol  spurious  magnetic  fields  in  the  balance  chamber. 

A  second  major  source  o£  spurious  recoil  force  was  inadequate 
shielding  of  the  "Photopot"  from  stray  light  in  general,  and  specifi¬ 
cally  irom  light  from  the  furnace.  Stray  light  from  the  furnace  i 

caused  a  shift  in  the  null  point  of  the  '’Photopot";  when  the  furnace 

power  was  terminated  to  measure  the  recoil  force,  the  stray  light  j 

from  the  furnace  was  also  terminated.  The  null  point  of  the  "Photopot" 
shifted  back  to  the  "furnace  cold"  position,  but  since  this  shift 

coincided  with  a  recoil  force  measurement,  the  shift  was  incorporated  j 

into  the  apparent  recoil  force.  Shielding  of  the  "Photopot"  was 
relatively  simple,  once  the  difficulty  had  been  identified. 
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The  residual  spurious  forces  which  remain  appear  to,  arise  from 
pressure-dependent  phenomena,  e.g.,  thermontolecular  flow  forces44. 
However,  these  residual  spurious  forces  now  constitute  only  —1$  of 
a  typical  expected  recoil  force,  compared  with  ~50  to  100$  before 
the  various  corrective  actions  were  taken. 

B.  EXPERIMENTAL  RESULTS 

We  have  only  recently  eliminated  (or  minimized)  the  spurious  recoil 
forces  discussed  above.  The  experimental  data  we  report  were  obtained 
with  the  modified  furnsce  and  automatic  control  system,  but  before  the 
spurious  forces  were  identified  and  eliminated,  end  are  presented  to 
indicate  the  precision  which  has  been  obtained  in  recoil  force  measurements 
with  the  microbalance. 

The  data  in  Table  7  for  the  vaporization  of  tin  were  obtained  with  a 
graphite  Miker  cell21  the  conical  orifice  of  which  is  described  by  T  - 

28.8°,  (L/r0)  »  4.98,  W  =  0.908.,  and  f  -  1.088. 

The  results  of  Table  7  arc  typical  illustrations  of  the  reproducibility 
obtainable  with  the  Miker  system.  It  should  be  noted  in  particular  that 
the  reproducibility  in  measurements  of  recoil  force  is  a_t  least  as  good 
as  the  reproducibility  in  measurements  of  rate  of  effusion. 

We  expect  to  observe  similar  reproducibility  end  good  agreement 
between  effusion  and  recoil  measurements  in  subsequent  vapor  pressure 
determinations. 
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TABLE  7 


HIKER  DATA  FOR  VAPORIZATION 
OF  TIN 


Tetap, 

Rate  of 
effusion, 
Ukit/aec 

Recoil 

Current, 

mA 

Recol 1 

Maaa*, 

.  HR. 

Iiawat, 

rE 

^  dro/cTr' 

PR 

1466 

1.061 

0.306 

100.4 

11.2 

21.7 

1470 

1.078 

0.310 

101.7 

11.4 

22.0 

1473 

1.138 

0.296 

97.7 

12.1 

21.1 

1473 

1.067 

0.302 

99.1 

11.4 

21.4 

1486 

1.156 

O.328 

107.6 

12.3 

23.2 

1489 

1.070 

0.327 

107.3 

11.4 

23.2 

Recoil  mast  la  that  tnaaa  which,  under  the  acceleration  of  gravity, 
counterbalances  the  recoil  force;  balance  calibration  ■  0.328  mg/mA. 


^The  "accepted"  value**  la  7.8  dyn/cra2  at  l480°K. 
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SECTION  VI 


CALORIMETRIC  STUDIES  OF 
VAPORIZATION  PROCESSES 

Determination  of  equilibrium  vapor  pressures  and/or  decomposition 
pressures,  with  subsequent  calculation  of  enthalpy  of  vaporization  and/or 
decomposition  and  of  enthalpy  of  formation  o£  the  gaseous  species  evolved, 
occupies  a  key  role  *  one  is  tempted  to  say  the  keystone  role  -  in  high 
temperature  chemistry  (HTC).  From  the  various  measurements  of  (supposedly) 
equilibrium  pressures  in  various  laboratories  have  arisen  numerous  dis¬ 
crepancies:  between  Knudsen  effusion  data  and  mass  spectroraetric  data; 
between  Knudsen  effusion  data  and  Langmuir  vaporization  data;  and  even 
between  two  sets  of  Knudsen  effusion  data,  both  of  which  were  obtained  by 
reputable  workers  using,  supposedly,  the  same  technique-  Attempts  to 
resolve  these  discrepancies  have  resulted  in,  among  other  things,  the 
formation  at  the  3966  Gordon  Research  Conference  on  1ITC,  of  a  committee 
to  study  establishment  of  vapor  pressure  standards  above  500° K  and  below 
10”3  atmosphere. 

In  addition  to  refinement  and  standardization  of  effusion  techniques, 
an  endeavor  to  which  the  previous  sections  of  this  report  is  devoted, 
new  approaches  are  needed.  A  more  or  less  obvious  one,  except  for 
experimental  difficulties  encountered  at  high  temperatures,  is  direct 
calorimetric  measurement  of  enthalpy  of  vaporization  and/or  decomposition. 
Sunner  and  MorawetZ46  (at  the  University  of  Lund,  Sweden)  ha,'e  studied 
intensively  the  problem  of  calorimetric  measurement:  of  heats  of  vaporization 
of  various  hydrocarbons  and  other  organic  compounds  at  25°C.  They  have 
successfully  measured  heats  of  vaporization  with  a  precision  and  accuracy 
better  than  0.1  kcal/mole  for  materials  with  vapor  pressures  as  low  as  10"4 
torr  at  25°C.  They  have  also  measured  differences  in  heat  of  vaporization 
as  a  function  of  the  effusion  orifice  geometry. 


In  conjunction  with  work  supported  by  another  contract  we  became 
interested  two  to  three  years  ago  In  the  possibility  of  direct  calorimetric 
measurement  of  enthalpies  of  decomposition,  e.g.  A(s)  —  +  C( .  Our 

initial  approach  to  the  direct  measurement  of  enthalpy  changes  was  essentially 
to  copy  the  furnace  and  calorimeter  configuration  of  Speros  and  Woodhouse *8 
'.'quantitative  differential  thermal  analysis"  (QDTA)  system47  and  to  make 
various  changes  in  the  electronic  control  circuity  and  the  output  sig>  al. 

In  evaluating  our  QDTA  system  ,  we  measured  enthalpies  of  fusion  of  several 
metals  with  a  precision  of  0.5-1. 5$  an<*  an  accuracy  of  in  the  range 

150-k50oC.  Subsequently  we  investigated  the  decorapostiion  of  PbC03,  ZnCO^, 
and  NHjCl;  with  QDTA  we  were  able  to  determine  with  reasonable  precision 
and  accuracy  the  enthalpy  change  for  these  decompositions,  including  the 
change  for  each  of  three  steps  in  the  decomposition  of  PbCO; 

Our  experience  with  the  QDTA  concept  hae  led  us  to  make  several  sig¬ 
nificant  modifications  in  the  furnace-calorimeter  configuration  and  in 
the  control-output  circuitry;  we  refer  to  the  new  system  as  "Differential 
Scanning  Calorimetry"  (DSC)40.  These  modifications  are  presently  (Summer, 
1967)  being  debugged  and  tested.  We  expect  that  with  this  new  DSC  w; 
shall  be  able  to  measure  in  the  range  ^0-1000°C,  endothermic  enthalpy 
changes  resulting  from  fusion,  vaporisation,  decomposition  and  phaaa 
transition,  and  perhaps  even  exothermic  enthalpies  of  reaction. 

It  is  our  further  opinion  thst  it  is  feasible,  within  the  present 
state-of-the-art,  to  design  and  build  a  DSC  system  which  would  provide 
direct  measurement  of  enthalpy  changes  in  the  range  1000-2000° K,  and 
perhaps  to  2500°K.  Such  a  system  would  make  available  to  "higher"  tem¬ 
perature  chemists  the  capabilities  presently  available  in  commercial 
equipment49  up  to  500°C  and  expected  in  our  DSC  system  up  to  1000°C. 

Another  virtue  of  the  DSC  technique  is  this!  with  the  one  assump¬ 
tion  that  the  rate  of  energy  input  to  the  calorimeter  is  proportional  to 
the  rate  of  the  (endothermic)  process  occurring  in  the  calorimeter,  the 
data  one  obtains  are  directly  interpretable  as  the  rate  of  the  process 
as  a  function  of  temperature  and  of  time.  It  is  then  rather  straight 
forward47' 50  to  obtain  rate  constants  ovtr  a  range  of  temperatures,  and 
hence  the  activation  energy,  from  a  single  l-to-J-hour  run.  With  the 
older  QDTA  system  we  have  obtained40  energies  of  activation  for  the  various 


steps  in  the  decomposition  of  PbC03,  and  for  the  homogeneous  decomposition 
(in  solution)  of  a  complex  organic  azo  compound. 

In  summary,  we  are  convinced  that  direct  calorimetric  measurement  of 
enthalpies  of  decomposition,  vaporization,  etc.,  are  presently  achievable 
to  1000’C,  and  with  a  relatively  modest  development  effort  could  be  achieved 
at  2000°C.  While  these  techniques  arc  not  likely  to  replace  effusion 
techniques  in  the  near  future,  they  can  provide  for  high  temperature 
processes  data  obtained  by  other-than-equilibrium  techniques;  such  supple¬ 
mentary  and  complementary  data  are  sorely  needed  in  many  high  temperature 
chemical  systems. 


TABLE  8 


IDENTIFICATION  OF  MAJOR  COMMERCIAL  COMPONENTS  OF  APPARATUS 


Item 

Amplifier,  Lock- in 
PAR  JB-5 

Controller,  C.A.T. 

Type  10877 

Lamp,  Projector 

8V,  SOW,  Type  13113C-04 
Philips  (Holland) 

Null  detector,  d.c. 

Type  9834-2 


Photopot 


Power  Supply, 

30V  d.c.,  QB28-1 

Pressure  Meter, 

Equibar  Type  120 

Pump,  diffusion,  oil, 

4  in.,  PMC- 720 

Pyrometer,  optical 

Photoelectric,  "Fhotomatic" 

Valve,  vacuum, 
variable  leak 
Cat.  No.  9101-M 

Silicontrol  Pulse  Unit 
Type  VS6332AF 


Trap,  liquid  nitrogen 

It  r> _  -  .;ll 

Vj  I  y  U“  bdiu 

Tuning  Fork 
Type  40 


Suppl ier 

Princeton  Applied  Research  Corp. 
Princeton  Junction,  N.J. 

Leeds  and  Northrup  Company 
4901  Stenton  Avenue 
Philadelphia  44,  Pcnna . 

Local  Photographic 
Supply  Shop 

Leeds  and  Northrup  Company 
4901  Stenton  Avenue 
Philadelphia  44,  Penna . 

Giannini  Controls  Corporation 
55  N.  Vernon  Avenue 
Pasadena,  California 

Sorensen  Products 
Raytheon  Co. 

S.  Norwalk,  Conn. 

Trans- Sonics,  Incorporated 
P.0.  Box  328 

Lexington  73,  Massachusetts 

Consolidated  Vacuum  Corporation 
1775  Mt.  Read  Blvd, 

Rochester  3,  New  York 

Pyrometer  Instrument  Co. 
Bergenfield,  N.J. 

Granville-: Phillips  Company 
5675  E,  Arapahoe  Avenue 
Boulder,  Colorado  80301 

VecTrol  Engineering  Div. 

Sprague  Electric  Co. 

P.O.  Box  1089 
Stanford,  Conn. 

Granville-Phillips  Co. 

5673  E.  Arapahoe  Ave. 

Boulder,  Colorado  80301 

American  Time  Products 
61-20  Woodaide  Ave. 

Woodside,  N.Y„  11377 
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